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ABSTRACT 


This thesis investigates the estimation of Lanchester 
attrition-rate coefficients in a deterministic aggregated 
combat model. Numerical values for the attrition-rate coef- 
ficients in the Lanchester-type model are taken to be the 
same as those for a corresponding continucus-time Markov- 
chain model, and maximum likelihood estimators are developed 
for these Markov-chain coefficients. A discussion of the 
historical background of Lanchester's equations preceeds the 
theoretical development of components of the attrition rate 
equations. Using the functional forms and procedures devel- 
oped by Gordon M. Clark in his doctoral thesis, a simple 
computer simulation program is developed for a short dura- 
tion battle consisting of homogeneous forces. Recommenda- 
tions are made +o modify this program to model heterogeneous 
forces in a Similar battle. A discussion of the theoretical 


background supports these recommendations. 
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I. INTRODUCTION 

The purpose cf this thesis is to investigate methods for 
estimating Lanchester attrition-rate coefficients, used in a 
deterministic aggregated combat model, from data generated 
by the detailed, high resolution Monte Carlo combat sinula- 
tion model called STAR (Simulation of Tactical Alternative 
Responses). An additional purpose is to construct a simple 
Lanchester-type combat model capable of using the estimated 
attrition rate coefficients to predict the outcome of STAR 
Simulations. 

STAR is a model which has been developed primarily by 
U.S. Army students attending the Naval Postgraduate School. 
It has recently gained acceptance in the Army Community as a 
Viable combat model, useful for studies such as the United 
States Army Mortar Study to be conducted in 1982. The cur- 
rent version of STAR is written in SIMSCRIPT and requires 
about 2 megabytes of core storage. The run time for one 
replication of STAR involving a company size force defending 
against a battalion size attacking force is approximately 
four CPU minutes on the IBM 3033 computer. 


The basic conceptual difficulties found in this thesis 


are Summarized in Figure 1. The numerical values for the 





attrition rate coefficients in the deterministic Lanchester- 
type model are assumed to be the same as those for a corre- 
sponding continuous-time Markov chain model. With this 
assumption, maximum Likelihood estimators (MLE's) can be 
developed for the Markov chain model coefficients. Specific 
data needed for the MLE computations is extracted From the 
high resolution Monte Carlo Simulation and used to compute 
"estimated" coefficient values. These values for the 
continuous +ime Markov chain model are then substituted into 
the "analogous" deterministic Lanchester-type model. 

Constructing a Lanchester-type attrition model such as 
that proposed in this thesis is of significant value. With 
respect to the hierarchy of models concept [{ Ref. 1], it 
could possibly fill the requirement as a ‘hook® or link 
between STAR and low resolution models that exist at higher 
levels in the pyramid cf models. However, possibly the most 
important contribution this model can make is as a filter 
model for STAR to reduce costs and time to run the many var- 
lations of a scenario in STAR. To explain further, consider 
the upcoming mortar study. To conduct such a study many 
runs/replications of a battle using the STAR model will have 
to be made. These numerous runs are required to compare 


various types of mortars, various types of force mixes, 
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This requires significant computer time at a significant 
cost. By using a Simpler analytical model which produces 
comparable output to a STAR run, many force mixes/weapon 
mixes/scenarios can be eliminated from further consideration 
at a minimal cost. STAR can then be used to concentrate 
only on the more promising options for force/weapons mixes. 

This thesis records the progress made in developing the 
small, simple, Lanchester-type combat model discussed above. 
It also lays the foundation for future work on this model. 
Chapter 2 of this thesis provides aa introduction to Lan- 
chester-type combat models. Chapter 3 outlines the develop- 
ment of the Deterministic Aggregated Model of STAR on the 
Apple computer (DAMSTAC), which represents the first step in 
building an analytical model comparable to STAR. Chapter 4 
explores various Lanchester attrition rate functional forms 
which may be used in the DAMSTAC or in future generation 
models. Chapter 5 details some desired enhancements of the 
model, while Chapter 6 discusses the work required to tran- 
sition from «he current homogeneous force DAMSTAC model to a 
more realistic, albeit more complicated heterogeneous model. 
Chapter 7 contains the concluding remarks. 

Appendix A lists the variables used in the computer pro- 


gram for the DAMSTAC model. Appendix B lists some notes of 
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interest concerning the program and nodel, which may not be 
apparent to the reader. Appendix C discusses the utility 
programs developed to support the model. Appendix D isa 
discussion of additional work done by the authors in the 
form of a computer model that explicitly plays the terrain 
effect. Appendix E consists of work done by two other stu- 
dents on a postprocessor which supports the computational 
requirements of the model and acts as an interface between 
STAR and the model. Finally, Appendix F contains the compu- 


ter listings of all programs written for this thesis. 
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II. HISTORICAL BACKGROUND 

Lanchester-type combat models are analytical models in 
which differential equations are used to express the rate of 
change of force levels to represent the effects of the 
attrition process. The concept of representing attrition 
effects in land warfare with such mathematical equations was 
originated by F. W. Lanchester in 1914 in order to guantita- 
tively justify concentration of forces. 

Lanchester hypothesized that under conditions of "modern 


warfare attrition between two homogeneous forces, each 


using "aimed" fire, could be modelled by: 


dxfdt = -ay (t) with x(O)= x, 


dy/dt ~bx (t) With y(O)= y 


9 


where x(t) and y(t) are the number of X and Y firers alive 
at time t, and a and b are positive constants that are 
called Lanchester attrition-rate coefficients. The constant 
‘at represents the effectiveness of an x firer in Killing Y 


targets per unit of time and 'b' represents the equivalent 
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effectiveness of the y firers. The equations written above 
are usually called "Lanchester's equations for modern war- 
fare". Clearly, these functional forms for attrition rates 
that Lanchester developed are directly proportional to the 
attrition rate coefficients and the number of firers from 
+he opposing force. The assumptions that must be made in 
order to use Lanchester's aimed fire equations to model cca- 
bat attrition are very restrictive. These assumptions can 
be stated as follows [Ref. 2]: 

1. Two homogeneous forces are engaged in combat; 

2. Each unit on either side iS within weapon range of all 

units on the other side; 

3. The effects of successive rounds in the target area 

are independent; 

YW. Each unit is sufficiently well aware of the location 

and condition of all enemy units so that it engages only 

live enemy units (one at a time) and kills them at a 

constant rate, which does not depend on the enemy force 

level. When an enemy target is killed, search begins for 

anew target, with the rate of acquiring a new enemy tar- 

get being independent of the enemy force level; 

5. Fire is uniformly distributed over surviving enemy 


waits. 
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Lanchester also hypothesized that when both sides used 


"area" fire, the attrition of forces could be modelled by: 


dx/dt ~axy 


dy/dt -bxy 
The first three conditions under which Lanchester's 
equations for area fire apply are identical to those for the 
equations for modern warfare, but assumptions 4 and 5 are 
replaced by { Ref. 3]: 
4. Each firing unit is aware only of the general area in 
which enemy forces are located and fires into this area 
Without feedback about the consequences ef its fire; 
5. Pire from surviving units is uniformly distributed 
over the area in which enemy forces are located; 
6. Each unit presents the same vulnerable area to enemy 
fire. This vulnerable area is much larger than the effec- 
tive area of a single round of enemy fire. Additionally, 
the number of hits required for a kill obeys a geometric 
probability law. 
The state equation that relates the force levels and 
does not ¢xplicitly contain time for the aimed-fire assump- 
tions is known as Lanchester's Square Law, 


b(x? =< X<) = aye - y?*). 
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The state equation for the area fire nodel yields 
Lanchester's Linear Law, 
bax, -Se) = aly, - yr 

From these state equations many important results can be 
deduced, for example, battle outcome prediction conditions. 
Further analysis of Lanchester's classical equations and 
their implications will not be done here. Suffice it to say 
that Lanchester's original work forms the basis for the ana- 
lytical combat model, with his simple model serving as a 
point of departure for the more complicated and elaborate 
analytical models discussed later in this thesis. Readers 
of this thesis and students who intend to continue with work 
initiated here should have a sound understanding of Lanches- 
ter's basic models before continuing further. 

Lanchester's equations for modern warfare and his equa- 
tions for area fire have, over the years, been thoroughly 
dissected and analyzed by students of dirtferential equations 
and students of analytical combat models. One such student 
was Gordon Clark, a doctoral student at Ohio State Univer- 
Sity in the late 1960's. He was well aware of the pvower 
that Simulation models such as Carmonette, a high resolution 
combat model, had for describing the complex process of con- 


bat. However, the major drawback to simulation, as he 
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perceived it, was that the method was slow and expensive. 
Clark proposed that the use of an additional model to inter- 
pret simulation results would greatly benefit Army studies 
which required numerous runs with the simulation model. de 
went on to develop this analytical model which he called the 
Combat Analysis Model (COMAN). This much Simpler and faster 
running model would greatly reduce computer execution time 
and costs, and could serve as a filter model for Carmonette. 
A similar approach to Clark's methodology in developing 
COMAN will be used here to develop the analytical model for 
STAR. The next chapter details the first steps taken to 


build the analytical combat model for STAR. 
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The aggregated force-on-force analytical model is an 
integral part of combat modelling and serves several 
significant purposes. With the high cost of computer time 
and the limited computer core storage available to most nil- 
itary users, the analytical combat model offers a viable 
alternative to detailed high resolution simulations. 

Analytical combat models can be used in a stand alone 
configuration or in conjunction with a high resolution sinu- 
lation. In the stand alone mode, the models are usually 
"transparent", or, in cther words, easy to understand. They 
are also small, quick but most importantly, they increase 
the user's understanding of the process of combat. In con- 
junction with high resclution combat simulations, the ana- 
lytical model offers, as Gordon Clark states in his doctoral 
thesis, the opportunity to "extrapolate the results of the 
combat simulation to predict the outcome from force nixes 
that were not explicitly simulated." In this way it can act 
as a filter to eliminate certain force mixes from further 
investigation. Also, Clark identifies another use for it as 
"an integral comvoonent of a large unit model to predict the 


outcome from individual unit actions." [ Ref. 4] 
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A major objective of this thesis was to build a simple 
aggregated combat model on the Apple II microcomputer. The 
input parameters for the model were estimated using the out- 
put data from a STAR run. This concept of estimating the 
attrition rate coefficients involved additional work in 
developing a data postprocessor compatible with STAR, and 
significant effort in the statistical analysis of the infor- 
mation provided by the postprocessor. 

There are several reasons why the Apple II computer was 
chosen for the development of the model rather than a main 
frame computer. First, the need exists for a ‘portable’ 
model that is easily accessible to the Army analyst. A desk 
top type system that can run a model quickly and with mini- 
mum setup time will be a valuable asset to the analyst. 

With the Apple II computer gaining wider acceptance in the 
Army for various uses, the accessibility of a model on such 
a system is greatly increased. 

Second, the nature of the model is such that it lends 
itsel= to a small system format. The aggregated form of the 
model reduces the required core storage from that required 
by a high resolution model. The concept of DAMSTAC is to 
serve as a filter model to weed out unacceptable force nixes 


or scenarios before turning to a high resolution model. The 
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idea here is to reduce the time and cost of running numerous 
iterations of a large, high resolution model on a main frame 
computer. 

Third, the cost of a computer system to Support the 
model is certainly economical. At a cost of under $3000, 
the Apple II computer system with a model such as DAMSTAC 
provides a powerful analytical tool for a minimal cost, cer- 
tainly in relation to the cost of a main frame computer. 

Fourth, the Apple II has the capability to provide addi- 
tional analysis functions on the output of the model. With 
such programs as Apple Plot and available commercial statis- 
tical and numerical analysis packages, most of the needed 
analytical work can be accomplished on the Apple II compu- 
ter. This can result in additional savings in time by 
eliminating the reguirement to transfer any data from one 
system to another. 

The laurels cast on the Apple II should not go without 
some discussion of its limitations. The Apple II is limited 
in capacity +o 64 kilcbytes of Random Access Memory storage 
Space. So, essentially the model must be able to work 
within the amount of core storage. There are methods of 
Working around this problem, if it arises, but this will 


usually create a new problem or limitation. For example, 
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some of the data or a part of the program can be kept on the 
storage diskette and put into computer memory only when it 
is needed. This saves space in the computer's memory Dut 
increases the run time of the model due to the time required 
to access the diskette to retrieve the required information. 
This raises another limitation of the Apple II--its run 
time. Compared to a main frame computer, the Apple II is 
very slow. This is, however, relatively speaking since con- 
puter time is normally measured in milliseconds and nanose- 
conds. For the analyst waiting for the output from DAMSTAC, 
the model on the Apple II would take about four to ten 
Minutes to run. This is probably acceptable for most 


analysts. 


A. SCENARIO 

The model portrays a battle between a RED tank battalion 
attacking a BLUE tank company in the 10 x 10 KM Fulda Gap 
(Keyhole) terrain. See Figure 2 for the terrain map. 

The BLUE defensive positions and the RED positions along 
with their attack routes (coordinating points) are provided 
in Table 1. The RED units attack along three preplanned 
attack routes toward the BLUE defensive positions. Grid 
coordinates for unit locations and coordinating points were 


computed as the centroid of the elements of the unit. [In 
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Star's data base, each tank has its own individual location. 

For example, BLUE unit 1 has four different locations asso- 

ciated with the it, each one corresponding to a tank in the 

unit. The centroid of the unit was taken to be the average 

of these grid coordinates. Thus, the initial location for 

BLUE unit 1 in Table 1 is the average of these coordinates: 
54150 100690 54150 100620 


54190 100520 54160 100790 


The rate of RED advance is a constant but specified by the 
user. BLUE units can move to alternate positions at a rate 
designated by the user. However, for this scenario, BLUE 
units remain stationary for the entire battle. 

To obtain the raw output data that would correctly 
reflect the above scenario, it was necessary to "turn off" 
many of the submodels of STAR. For example, the Field 
Artillery module was shut off by scheduling the first fire 
mission well after the battle ended. Weapon systems other 
than tanks within the fighting units were rendered useless 
by setting their ammunition counts (basic loads) to zero. 
Additionally, the communications and the Air Defense modules 
were not played. After all of these preliminary actions 
were taken, the result was a battle exclusively between 


tanks. This effort worked well to reduce a complicated, 


24 





heterogeneous force battle to a simple homogeneous force 
battle. One minor problem arose, however, when a weapon 
system other than a tank was detected and fired on when, by 
+he scenario, it did not even exist. The effects of this 
target acquisition and engagement were considered insignifi- 


cant and were disregarded. 


Table 1: UNIT LOCATIONS/COORDINATING POINTS 


BLUE X-COORD Y-COORD SIZE REMARKS 

UNIT 1 54163 100655 i initial location 

UNIT 2 54788 98120 ~ mnatwal location 

UNIT 3 57448 97893 4 Pitt ale Ocatton 

RED X= COORD Y=COORD SIZE REMARKS 

UNIT 1 98727 101491 70 Llnitial location 
52000 96000 coordinating pt 

UNIT 2 58110 10 16 96 10 Special Locatson 
52000 96000 coordinating pt 

UNIT 3 58727 101491 10 tata di LocatTon 
52000 96000 coordinating pt 


B. ASSUMPTIONS 
The attrition rate functional form chosen for the 
DAMSTAC model, which will be discussed later in this thesis, 
requires the following assumptions. 
1. The conditional kill rates, @~and B , are constant and 
independent of time if an individual firer has at least 


On? acquired target. If there are no acquired targets, 


Zo 





the kill rates are zero. 

2. Pandq, the probabilities of nondetection, are 

assumed to be constant functions of time, that they are 

independent of the force size and the number of acquired 

targets, and that the probabilities p and gq are equal for 

all x-firer/y-target and y-firer/x-target combinations 

respectively. 

3. Synergistic effects are neglected, therefore, the 

actual force attrition rate is equal to the sum of the 

individual kill rates for the opposing forces. 
Additionally, the assumptions listed below are designed to 
remove the distractions of a complicated scenario and to 
make the initial model easier to understand: 

1. All units have point locations; ie. a battalion is 

represented by a grid coordinate and all calculations for 

the battalion are made from that location. 

2. Units are homogeneous. That is, they are composed of 

identical firing elements/weapons systems. 

3. Fire is uniformly distributed over surviving enemy tar- 

gets within range. 

4. Sufficient supplies of ammunition for both sides exist 

mec the duration of the battle. 


5. Units move at a constant rate of speed. 
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C. SYNOPSIS OF THE COMPUTER MODEL 

The DAMSTAC computer program is a modular type program 
writcen in the Applesoft BASIC language. The program 
accesses external data files as a source of input of initial 
data for the two forces. These data files are produced by 
the user with the help of the Force Maker utility program 
described in Appendix C. Additionally, the program has the 
capability to make output files that can be stored on a 
diskette and accessed by the user with the Results Reader 
utility program, also described in Appendix C. Results 
Reader will print the output to the screen or +o a printer. 

A copy of the program for the model is included as List- 
ing A in Appendix F. A sample run of the model including 
user inputs is included as Figure 3. Flow charts of the 
main program and subroutines are included as Figure 4, Fig- 
ure 6, and Figure 7. 


1. Initialization of Arrays and Variabies 





Lines 30 thru 170 dimension the arrays and matrices 
and initialize hardwired data. Dimensions are set to handle 
up to five units in each force and up to ten coordinating 
points, including the initial location, for each unit. The 
user can modify these capacities by changing the dimensions 


of the appropriate variables. 
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Lines 180 thru 350 read BLUE force data from a file 
to include the number of units, initial locations, any coor- 
dinating points, and rates of advance, and compute total 
force level based on user input of the unit force levels. 
The time increment (DT) is set to 30 seconds while the maxi- 
mum time of the battle (MT) is set to 1500 seconds. These 
values are, of course, adjustable by changing them in line 
80. For both forces, the breakpoint criteria is set at 0.3 
of the force level and movement criteria is set at 0.5 of 
the unit strength, even though the movement criteria was not 
used for any runs of the model. 

Lines 360 thru 520 perform the same operations for 
the RED forces as lines 180 thru 350 perform for the BLUE 
forces. 

renee 530 and 600 ask if the printer should be 
turned on or off and calls the printout subroutine to print 
the initial data. 

2. Main Program (lines 610-780) 
The main program performs these operations: 
increments the time, 
checks if the time exceeds the max time, 
calls the movement and attrition subroutines, 


recomputes current total force levels after attrition, and 
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checks force levels against the designated breakpoints for 
zach forcs, to determine it the pattle has ended. 
If the breakpoints or the max apy ARE exceeded, the battle 
terminates and the program prints the force levels. Addi- 
tionally «he program will write the force levels for each 
sime interval onto the diskette. (lines 1470 thru 1580) 


3. Movemen* Subroutine (Lines 790-1160) 





Both RED and BLUE forces are capable of moving. 
Movement of the BLUE force is nullified by setting its nove- 
men*= rate to zero in the initialization part of the progran. 
The rovement subroutine uses basic geometry to compute new 
coordinates for «he moving unit (see Pigure 5). The nove- 
ment process is based cn noving a unit from one coordinating 
point co the next in a straight line. The amount of nove- 


ment depends on the rate of movement (BR(I) or RR&R(I)) and 


NEXT 
CaoRO. PT. 


CURRENT 
UNIT LOC. 





Figure 5: Geometric Representation of Unit Movement 
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(SEE FIGURE 5) 
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AS DIST. To 
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USE MOVE 






COMPUTE 
NEW 





LOCATION 


Figure 5: Flow Chart of the Movement Subroutine 
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the time interval length (DT). If the distance computed is 
farther than the next coordinating point, only movement to 
the coordinating point is accomplished. Lines 800 to 970 
perform computations to move RED forces while Lines 980 to 
1150 perform the same computations for BLUE fcrces. 


a. ttrition Subroutine (lines 1170-1480) 





A range check is made for each firer/target combina- 
tion to insure the units are within firing range. [If so, 
attrition caused by that firer on the target is computed. 
Next total attrition for each unit is computed based on the 
partial attritions of the firers as a fraction of their 
force levels. In other words, when a firer is firing at 3 
enemy units, each target unit receives 1/3 of the firer's 


fire power. 


There are two different output formats used in the 
program. The first (lines 1490 to 1530) prints only a sun- 
mary of the total force leveis for both forces for each time 
interval. This is the same data that is saved to a file by 
the final portion of the program. The second format (lines 


1760 +0 1940) prints out the units of seach force with their 
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Pigur2 7: Flow Chart of the Attrition Subroutine 
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current locations and strengths. This format allows the 
user to see which unit, specifically, is being attritted and 
the grid coordinates of that unit at the time of attrition. 


6. OQutput File Making Subroutine (Lines 1540-1750) 





This subroutine allows the user to make a file of 
the output from the current run of the program. This 
routine actually produces two files, one consisting of the 
attrition rate coefficients and probabilities of non-acqui- 
Sition, and the second containing the time and force levels 
of both forces at that time. These files are seperated so 
that the second file can be used with the APPLEPLOT program 
to produce a graph of the two force levels as a function of 
time. This file has the postscript, ".RESULTS", added to 
its file name, while the first file with the attrition rate 


coefficients has the postscript ".COEFS". 


D. METHODOLOGY FOR COMPUTING PARAMETER ESTIMATES 

The process of building a computer simulation combat 
model requires the deveiopment of methodologies for 
modelling the various aspects of combat. Constructing a 
force-on-force attrition model requires methods for consoli- 
dating combat phenomena into simplified algorithms with 
coefficients representing one or more facets of combat. 


This section discusses one significant part of the 
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force-on-force attrition model: the methodology for estimat- 
ing the attrition rate coefficients in the Lanchester type 
attrition rate equations. 

In deriving *the equations for the attrition rates (dx/dt 
and dy/dt), it is necessary to understand the processes that 
must occur to cause a firer to inflict attrition. These key 
processes can be summarized as: 

the environmental process; 

the target acquisition process; 

the target engagement process. 
The environmental process involves the terrain and weather 
factors that affect firer/target line of sight (LOS). The 
target acquisition process represents the visible capability 
of the firer to detect a target on the battlefield, while 
the target engagement process describes the firer's ability 
+o hit and kill an acguired target. These three processes 
will all be represented by the attrition rate coefficient, 
A, in the attrition rate equation. This coefficient will be 
broken down into its three components which describe mathe- 
Matically the phenomena of the three processes. The three 
components are: 

Pp = the probability that a target is visible to the firer 


V 


(environmental process) 
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P, = the probability that a visible target is acquired 
(target acquisition process) 
@= the casualty rate at which a firer attrits the 
acquired target (target engagement process) 

a should not be confused with A, the attrition rate 
coefficient, of which q is a factor. The attrition rate 
equation may be represented in words by: 
dx/fdt = -[P(tgt is visible)][P(tgt is acquired)][casualty rate] 
The second probability above ae be broken down further to 
P(tgt is acquired | tgt is visible). ‘Note that this is a 
conditional probability since it is assumed that the target 
must be visible to be acguired. By identifying the 
mathematical expressions for each of these components of the 
attrition rate coefficient, the attrition rate equation can 
be obtained. 

1. The Environmental Process 

There are generally two ways to represent line of 
Sight in a model: 
as a mathematical simulation of actual terrain; 
as a stochastic process. 
The mathematical simulation method normally involves a 


digital representation of terrain elevation over a segment 


of the terrain. The elevations of the segments between the 
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firer and the target are used to determine if line of sight 
exists. Virtually all high resolution models use this 
method, or variations of it, to determine line of sight. 
STAR uses a different form of mathematical simulation to 
represent terrain. This is called functional terrain, in 
which hills and forests are represented by elliptical con- 
tours of elevation. Hills are modelled by describing then 
in terms of bivariate normal equations with varying parame- 
ters for each hill. Forests or built-up areas are repre- 
sented by an elliptical pattern on the ground again with 
varying parameters for each forest. 

For a low resoluticn or highly aggregated model, the 
mathematical simulation method for line of sight is usually 
not practical or desireable in terms of core space and con- 
putation time required. Therefore, the second method, in 
which LOS is represented as a stochastic process, is used. 
The concept of this process is as follows. The effect of 
terrain on the process of line of sight is that a target 
falis into one of two categories or states: either it is 
visible or it is invisible. 

Assuming a steady-state situation for the terrain, 
probabilities can be assigned to each of these two states. 


These probabilities represent the percentage of time that 
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the target will be in that state (visible or invisible) for 
that terrain over a reasonable amount of time. Graphically, 
the states of the target can be represented by the diagran 
fae fagucre 8. 


TARGET 
Visit @Lée 


TARGET 
tWVisi Be 





8 
TIME Ts 


Figure 8: Two-statse Stochastic Graph of Intervisibility 


Since this is a stochastic process, there is an 
exponential distribution associated with each of the two 
states. This distribution represents the time between two 
occurtences of the same state. On the graph in Figure 8, 
*he time between the end of rm and the start of oa 2s expo- 
nentiaily distributed with mean, 7. Likewise, the «ine bet- 
ween *he end of te and the start of Tis exponentially 
distributed with mean, uw. Graphically, this can be repre- 
sented by Figure 9, where the circles tepresent the «wo 
states, visibility and invisibility, and the arrows 
Tepresent “he transition from one state to the other at the 


given rats. 
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VISIBLE INVISIBLE 
STATE STATE 


Figure 9: Two-state Stochastic Model of Intervisibility 


The attrition caused by a firer on a target can 
occur only when the target is visible (assuming no attrition 
from wild, blind shots or indirect fire). Therefore, «he 
attrition rate is a function of the probability that the 
target is visible *o the firer. This probability will be 
denoted as Poe and can be expressed as the proportion of 
time that the target is in the visible state in a steady 


state situation, LS 


7 + p 


2. The Target Acquisition Process 
Assuming now that line of sight exists between a 

Sizer and a target, the process of target acquisition must 
now 3e@ represented. There are two primary modes of target 
acquisition: 

Serial acquisition, and 

PaedilLel acquisition. 
Tn the serial acquisition process, targets are acquired 


between engagements. That is, in concept, the firer is 
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concantrating on the engaged target and does not or is not 
able +o acquire other targets on the battlefield. Any addi- 
tional acquisitions made prior to engaging the current tar- 
get are lost and must be re-acquired after the current 
engagement has ended. Examples of firers who use serial 
acquisition are; 

any firer who is "buttoned up" and has his field of view 

narrowed to a thin band; 

wire guided missiles/anti-tank systems. 

In the parallel acquisition mode, the firer main- 
tains his list of acquired targets throughout an engagement 
(except those lost through loss of line of sight). Thus, 
when an engagement has ended, by the target being killed or 
loss of line of sight, the firer may immediately engage a 
target that was previously acquired and is still a valid 
target. All systems which do not use Serial acquisition are 
assumed to use the parallel mode. 

For purposes of this thesis, all firers are assumed 
to use parallel acquisition. Therefore, serial acquisition 
will not be discussed any further. The reader should note 
that the serial acquisition process does have validity in 
modelling action on the battlefield, and should be consid- 


ered in future enhancements of this model. 
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To explain the derivation of the target acquisition 
portion of the attrition rate equation, the diagram in Fig- 
ure 9 must be modified to chat in Figure 10. 


VISIBLE AND 
ACQUIRED 


STATE 
A 








INVISIBLE 
STATE 


VISIBLE AND 
NOT 


ACQUIRED 
STATE 


Figure 10: Three-state Stochastic Model of Intervisibility 


The «wo-state stochastic process has been changed to 
a three-state process by dividing the visible state into the 
visible-not acquired state (VNA) and the visible-acquired 
state (VA). What has cccurred is that a target can move 
from the visible-unacguired state to the visible acquired 
State at a rate A. Since both of these states comprise che 
Visibility state in total, the rate of moving from either of 
them to the invisibility state remains atu. Notice that 
there is no direct «transition from the invisible stata to 
“he visible-acquired state. This is only a matter of 
representation, since it is assumed that no more than one 


event, ie. moving to a new state, will occur during a tine 
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increment. Thus, the situation described can still be 
represented by Figure 10. 

Another point of interest in this diagram is the 
fact that there is no transition from visible-acquired to 
visible-not acquired. This implies the assumption that once 
a target is acquired, it remains acquired unless it is 
killed or line of Sight is lost, in which case it transfers 
to the invisible state. 

To determine the effect of the target acquisition 
process on the the attrition rate, the probability that a 
target is acquired must be determined. From Figure 10, and 
uSing the same concept of a stochastic process as that used 
in section III-D-1, the equation can be written 

P_(tgt acquired) = 4 : (1) 
Thie should more properly be stated as the probability that 
a target is acquired given that it is Visible, since this is 
actually a conditional probability with the state of acqui- 
Sition being dependent on the target being visible. 


Therefore, the probability that a target is visible 


and acquired in the steady state is: 


rg 
il 


a P(tgt vis & acq ) 


P, (tgt acgq |tgt vis) Pi (tgt vis) (2) 


Gar) rt) 





The implications of this equation will be discussed later in 
Chapter V. 
3. Formulating the Attrition-rate Equation 
Using equation (2) in the attrition-rate equation, 


the attrition rate of a single firer on a single target can 


be summarized using expected values as: 


dx/fdt = -P a (3) 


That is, the attrition rate equals the rate of killing an 
acquired target times the probability that a target is visi- 
ble and acquired. However, Since there are likely to be 
more than one firer and more than one target on the battle- 
field, this equation must be modified to reflect this. 

If a firer is going to cause attrition on his tar- 
gets, he must be able to acquire one or more of these tar- 
gets. For firer yl with x targets, the probability that y1 
acquires one or more targets is just one minus the probabil- 
ity that he acquires no targets. Assuming the same 
probability of acquiring any of the targets, the probability 
x 


of acquiring none of the targets is (1 - P 


a - Therefore, 


P = > iG Po) (4) 
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The attrition rate now becomes 


dx/dt = -P a (5) 


for a single firer against x targets. Since there are more 
than one firer on the battlefield, each assumed to be acting 
independent in this process, the attrition rate in equation 
(5) must be multiplied by the number of firers(y) and the 
equation is now complete; 


dxsdt = -P aye (6) 


In terms of the probability of acquiring a target, this 


equation can be written as: 


iyde == {1 = Pe) ) ay. (7) 


Va 
This is the same equation which Clark uses in his COMAN 
model, that uses results from a high resolution model to 
estimate the value for the conditional casualty rate and the 
probability that a target is not acquired. Henceforth, (1 - 
Fae Will be referred to as p for RED firer/BLUE target con- 
binations (or q for a BLUE firer/RED target combination). 


lark) 


3 


he target engagement process can be modelled in its 


most detailed form, as in the high resolution model, or in 
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its simplest form, as in the aggregated, analytical model. 
The high resolution model, such as STAR, includes explicit 
representation of every factor that is involved in the pro- 
cess of aiming at, firing at, hitting, and killing a target. 
These factors include: 

weapon system characteristics; 

weapon/fammunition type characteristics; 

target characteristics such as size, shape, vulnerability, 

etc. 

In contrast, the low resolution, aggregated model 
represents the engagement process by one value--the attri- 
tion rate coefficient. This is the method to be employed in 
the model presented here. The task now is to develop a 
method to accurately estimate this coefficient for use in 
the attrition rate equation. Note that we are actually 
talking about two attrition rate coefficients and two equa- 
tions. The derivation and method for eStimating one is the 
Same as the other, so for purposes of the discussion, refer- 
ences to an attrition rate coefficient will be in the 
Singular. 

There are presently two basic approaches to the 
problem of estimating attrition rate coefficients for paral- 


lel acquisition and stochastic line of sight. The first is 
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the use of the independent analytical model, involving 
weapon system characteristics and the use of algorithms to 
determine the expected time for a firer to kill a single 
target. The attrition coefficients are then the reciprocal 
of these expected times. This method will be referred to as 
Bonder's method after Seth Bonder who has used it exten- 
sively in his combat models such as VECTOR II. 

The second method is the fitted-parameter analytical 
model, involving the use of a separate high resolution, 
Monte Carlo combat simulation to produce data to determine a 
maximum likelihood estimate of the attrition coefficients. 
This method will be referred to as Clark's method, since 
Gordon Clark popularized it in his doctoral thesis and in 
the development of the Combat Analysis Model (COMAN). 

Most of this section involves a discussion of 
Clark's method and the practical application of it. How- 
ever, in order to provide a contrast to that method, a short 
description of Bonder's method will be provided. 

In general, the concept of Bonder's method is to 
produce an algoritha representing the expected time that it 
takes a firer to kill a single target. The elements of this 
algorithm include computations for line of sight, target 


acquisition, target selection, and the killing process. 
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Bach element accounts for a certain amount of time in the 
expected time to kill a target. This entire process is seen 
by Bonder as a Markov renewal process where the outcome of a 
round fired is dependent only on the outcome of the previous 
round. By using weapon systems characteristics data as raw 
input, times to acquire, select, and kill a target can be 
determined and used in the algorithm for expected time to 
kill a target. 

Bonder*s method uses an independent, analytical 
model requiring a large data base of factors, such as range 
dependent weapon system characteristics, gquantifications of 
line of sight and target acquisition parameters, to support 
the computation of the attrition rate coefficient. 

A detailed explanation of Bonder's method for esti- 
Mating the attrition rate coefficients can be found in the 
basic documentation for the VECTORII model [Ref. 5]. Basi- 
cally, Bonder says that the attrition rate coefficient A,.is 


J 
represented in the parallel acquisition mode by Pe On y 


B 1) 
where {> the probability that a given group-i weapon is 
firing at a group-j target. Q. 5 can be expressed in terms 
of the probability that a group-j target is visible and has 


been detected by a group-i firer. This probabiltiy is 


dénoted as Sia and the expression for Q.. is: 
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Q:.5 sti, - $,;) for 422. 
Assuming homogeneous forces, this equation reduces to 

Q;. = ane peanGe: ne 2— 1) Subscripts can be eliminated 
since there is no distinction among types of force elements. 
His resulting expression for S is identical to that of equa- 
tion 4 , where S35 and Sa ee defined the same. 

Bonder's expression for the conditional casuaity 
rate, @ , iS written as the inverse of the Pesckes time to 
kill an acquired target. This expression, when further bro- 
ken down contains elements dependent On the weapon system 
characteristics, target characteristics, and a number of 
other subsystem performance parameters. These parameters 
were determined empirically outside the operational environ- 
ment of the battlefield and molded together in this pre- 
determined structure to produce the value of the attrition 
coefficient. This derivation of the attrition coefficient 
is the point where Bonder's method and Clark's method 
ieerer. 

Clark's methodology is based on the assumption that 
the process of attrition in combat is a nonhomogeneous, two 
independent type Poisson Process where the time between 


casualties has an exponential distribution with parametera”, 


Called the total conditional casualty rate. Likewise, the 
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time between RED casualties and the time between BLUE 
casualties each has an exponential distribution with parame- 
ters A+ and , respectively, where A =A,tApL- From the 
theory of Poisson Processes, the probability of a RED casu- 
alty occurring isAr/, , and the probability of a BLUE casu- 
alty occurring is 5/,, - With this basis, Clark is able to 
develop a likelihood functon in terms of the conditional 
casualty rates q@ and gp , of which the casualty rates are 
functions, ie., Rime f(a) and Ay= £(¢). (Ref. 6] It should 
be noted that, in order to maintain consistency of terminol- 
ogy, the equations attributable to Clark's work have been 
adjusted to reflect the definitions presented in this the- 
Sis. 

Clark also uses the concept that attrition is a 
function of target acguisition, which depends on the prob- 
ability of a firer acquiring a target. This concept pro- 
duces two more parameters which must be estimated: p, the 
probability that a BLUE target is not acquired by a RED 
firer; q, the probability that a RED target is not acquired 
by a BLUE firer. 

For the functional form in Clark's COMAN model, the 


equations can be written as: 
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f= -afi-p™] with m(0) = M, 
- (8) 
fp= -B{i-q°] With no)= no 


where the parameters are defined by: 


f, = the conditional attrition rate function for BLUE 
forces 
f= the conditional attrition rate function for RED forces 
@= conditional kill rate for a RED weapon against acquired 
BLUE targets 


B= BLUE weapon conditional kill rate corresponding toa 


p = probability that a BLUE target is unacquired by an 


individual RED firer 


q = probability that a RED target is unacquired by an 


individual BLUE firer 
m = number of BLUE firers surviving 
nh = number of RED firers surviving 


These functional forms depend directly on the abii- 
ity of firers to detect targets. To see this, consider just 
Mmemconditional attrition rate function for the BLUE force, 
f,- The definition of p implies that the value of (1-p™) is 


the probability that a RED firer has at least one: acquired 
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BLUE «arget. Therefore, the attrition rate for the BLUE 
force is the kill rate for acquired targets times the prob- 
ability of having an acquired target times the number of 
firers. 

It should be noted that Lanchester's classical 
attrition rate functional forms are special cases of Clark's 
attrition rate functional forms, and, under appropriate con- 
ditions, Clark's equations reduce to Lanchester's equations. 


When p and q are zero, or very nearly zero, 


- @ (1-0O")n =- an (9) 


th 
> 
it 


and, 


- 6(1-0")m = - Bo (10) 


rh 
4 
i! 


which are Lanchester's classical aimed fire equations. When 
individual targets become increasingly difficult to acquire, 
p and gq assume values close to one and it can be proven that 


in the limit: 
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ae = - a@ (1-p) an (11) 


and, 


f. = - 8 (1-q)amn (12) 


Meetme attrition rate coefficient ais set equal to a({i - 
a end the coefficient b is set equal to g(1- q), then the 
two equations above are identical to Lanchester's equations 
for area fire. 

In his doctoral dissertation, Clark developed the 
likelihood function for a heterogeneous nodel, however since 
“he DAMSTAC model assumes homogeneous forces, the necessary 
equations have been reduced to a siaupler form than that 
derived by Clark. By differentiating the logarithm of the 
likelihood function with respect to p, the following equa- 
sion is derived: 


Saen ( >. @) oa Cama DoT L 


= eer CL - eal By 
al tate t%) 29 03) 
where 
C. = 0 if casualty k is RED 
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1 if casualty k is BLUE 


ae number of BLUE survivors pricr to kth casuait 
ay number of RED survivors prior to kth casualty 
ae eee, = time between casualty k & casualty k-1 


Likewise, by differentiating the logarithm of the Likelihood 
function with respect to @, and by solving for @ , a second 


equation in terms of $ and @ is produced: 


‘ a (14) 
(6 4 


1 a i me CN ra oa ae ) 


Differentiating with respect to g and B produces equivalent 
“ype equations which will allow the estimates for g and B to 
be determined in the same on as follows. 

Generally, the two equations (13) and (14) can now 
be solved Simultaneously to determine the values for 5 
and @. Uneeprinacsly, this 1s ¢aSier said than done. The 
complexity of these equations (more so in the hneterogeneous 
case) does not allow explicit equations for 5D and @. The 
semainder or this chapter will discuss. a nethod for deter- 
Mining the astimates of p and @ from these equations and 2 


way [Oo @xtract the necessary data from a nigh resolution 


Monts Carlo simulation run to compute these estimates. 
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E. ESTIMATING PARAMETERS FROM A HIGH RESOLUTION MODEL 
As mentioned before, equations (13) and (14) must be 
solved simultaneously to produce values for @ and D. Gener- 
ally, these equations are too complicated to produce such 
values easily so other techniques must be used. The method 
suggested by Clark is to enter a value for p in equation 
(14) and solve for @. This value and the selected value 
for 35 are then put in equation (13) and the equation is 
evaluated. yeuves that both values are estimates since we 
are computing the maximum likelihood estimates of p and gq. 
The values of p and @ which give a value for equation (13) 
closest to zero are the best maximum likelihood estimates 
for p and @~@. If equation (13) equalled zero exactly, the 
true values of the maximum likelihood estimates would be 
known. 
In uSing these equations, the input received from the 

high resolution, Monte Carlo Simulation gust include: 

the time between casualties; 

the number of RED and BLUE survivors prior to the kth 

Sasualty; 

the total number of casualties for RED and BLUE for the 

battle; 


the type of casualty (RED or BLUE) of the kth casualty. 
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In Clark's theory for a homogeneous force battle, the 
conditional attrition rates and the probabilities that a 
+arget is not acquired are assumed to be constant over the 
entire battle. This means that equations (13) and (14) are 
summed over all casualties in the battle and the estimates 
for p and a, which are determined from them, are also 
assumed to be constant over the entire battle. 

In applying this theory to STAR output, the authors 
obtained guestionable results for the values of the parame- 
ter estimates. The estimates, D and a were computed to 
have values of 0.99 and 0.93 respectively. These values 
imply that both BLUE and RED forces experienced difficulty 
in acquiring targets. Analysis of the Fulda Gap terrain and 
consideration of the respective locations of the opposing 
forces support the estimated values for p and g. However, 
the stated value for De 0.99, 1S a truncation of the actual 
value computed. If proper rounding procedures had been 
used, 5 would have been set equal to 1.0, which means that 
the RED forces would not detect any BLUE targets. Nondetec- 
*ion implies that no BLUE attrition will occur. But this 
Semtradicts the STAR battle which did, in fact, produce BLUE 
casualties amounting to about two-thirds of the total BLUE 


POrce. 





The estimates, aand B were computed to be 2.4 and 0.03, 
respectively. The relatively large value of @ implies that 
when the RED forces are able to detect BLUE targets, the 
BLUE eyecacaiinis suffer heavy casualties. There appears to 
be an unusually large difference between the values for da 
and B which are not reflected by the values for § and qe 
This disparity is most likely attributable to the restric- 
tion that all parameter values be held constant over the 
entire battle. Experimentation with different parameter 
values in DAMSTAC indicates that a value of about 0.40 more 
closely approximates the attrition observed in the STAR bat- 
tle. Experience with these computations reveals that this 
is not necessarily the case and that a more likely assump- 
tion is that the parameters are constant over different time 
intervals of the battle. In general, a simple battle can be 
divided into phases where the intensity of the battle is 
constant within each phase but different among the phases. 
Basically, the phases are: 

movement to contact (low intensity) 

initial contact (mid intensity) 

Main battle (high intensity) 

withdrawal from contact or annihilation (mid/low inten- 


sity) 
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Moving from one phase of the battle to the next should cause 
a change in conditional casualty rates (qaandf) as well as 
+he probability of not acquiring a target (p and q). 

The question now becomes how to determine the breakdown 
of the battle by phases using the output from a high resolu- 
tion, Monte Carlo simulation. That is, where do the changes 
in phases occur in the simulation? One possible method is 
to look at the number of casualties as a function of time. 
Another is to look at the number of shots fired as a func- 
tion of time (see Figure 11). In both cases, time is seg- 
mentei into time intervals (eg. every 100 seconds) anda 
histogram of the casualties or shots is plotted for each of 
these intervals. Having done this for several battles, the 
results identify fairly clear choices for the end of one 
phase of the battle and the beginning of the next phase. 
However, mor analysis of these two methods reveals that 
*hey night not truly identify the phases of the battle and 
shat an alternate method should be considered, primarily one 
based on range. This concept is discussed in detail in 
Chapter VI. 

Using the number of casualties per time interval, a dis- 


tortel picture of the intensity of the battle might be 


obtained. The interaction of the probability of acquiring a 
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Pigur2 11: Histogram of Shots per Time Interval 


Buege. With «he at«ricion tate of a fizrer might offset each 
other and produce results which conflict with the true bac- 
*le iatensity. In cther words, if the acquisition probabil- 
ity is high but the attrition rate is low, the number of 
casualties in a time interval night indicate a low intensity 
Dhase of «he battle when, in fact, it is a high intensity 
Phase. Therefore, this method should be used with caution, 
if at all. The method of using the number of shots per tine 
interval has similar faults. Twenty shots by one firer in a 
time interval is much different in terms of intensity than 
one shot by twenty firers. Again, this method should be 


used with caution. 
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The question is then: what method should be used? One 
method, which is used by the Combined Arms Combat Develop- 
ment Activity (CACDA/CASAA), is meaSuring the number of 
shots per firer per time interval.! Although it also has 
some faults, it does provide a better indication of the 
intensity of the battle. It is, for the most part, indepen- 
dent of the conditional casualty rates, although the number 
of shots would indicate a dependency on target acquisition. 
This dependency is reduced by averaging the number of shots 
over the number of firers. In effect, this measure can be 
considered to be unbiased by the changes in conditional 
casualty rates and probabilities of not acquiring targets. 
Thus one is able to get a feel for the intensity of the bat- 
tle from the histogram of the number of shots per firer per 
ime interval. This 1S not necessarily scientific reason- 
ing, but it seems to work. 

Once the phases of the battle are delineated for a high 
resolution, Monte Carle simulation run, the next step is to 
determine how to treat the 'dead* time from the beginning of 
a phase to the first casuaity and from the time of the last 
casualty to the end of the phase. For phases in the middle 


isonversation with Mr. Ernest Boehner from CASAA on 5 
December 1981 
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of the battle, each of these two times is a portion of the 
time between two casualties and should thus have a bearing 
on our maximum likelihood estimates from equations (13) and 
(14). 

To resolve this situation, the time of the last casuaity 
within a phase is used as the end of the phase. For exanm- 
ple, if a phase is supposed to end at 1200 seconds and casu- 
alties occur at 1140 and 1225 seconds, the phase will be 
+erminated at 1140 and the new phase begun at 1141. This 
also allows this time between casualties to be accounted for 
in the next phase. This method produces some variation in 
phase lengths when running several replications of the bat- 
tle, but this variation is minimal, especially since the 
beginning and end of phases are being judged from something 
less than a scientific method. 

S2veral replications, as many as affordable, should be 
used to compute the maximum likelihood estimates for the 
four parameters. By using a large sample population {the 
Casualties), a more accurate eStimate of the parameter can 
be obtained. Likewise, when the heterogeneous model is 
used, a larger sample population is desirable, since the 
number of engagements between two particular weapons systems 


will be much less than the number of engagements in the 
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homogeneous case. This could produce values that are 
statistically insignificant. ; 

A postprocessor to the high resolution model was devel- 
oped to make the numerous calculations required to produce 
these estimators. The postprocessor uses an optimization 
routine to determine the maximum likelihood estimates fron 


the given equations. Detailed information on the postpro- 


cessor is provided in Appendix E. 
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Different military situations have been hypothesized to 
yield different functional forms for Lanchester-type equa- 
tions. The term "Lanchester-type equations" refers to any 
differential-equation model of combat and the tern "func- 
tional form" relates to the structure of the attrition rate, 
i.e., how the attrition rate equation is written and which 
parameters are considered explicitly in the equation. 
Henceforth, a standard shorthand notation will be adopted to 
refer to the functional forms used in homogeneous-force Lan- 
chester-type combat models. For example, Lanchester's equa- 
tions of modern warfare are expressed as FIF functional 
forms (or, simply FIP attrition) because they are propor- 
tional to only the number of enemy firers. Taylor [Ref. 7] 
provides a good explanation of this convenient shorthand 
notation. 

Generally, the classical Lanchester functional forms for 
expressing the attrition rate for forces using aimed fire or 
area fire are tco simple to be used "as is" ina detaiied 
analytical combat model. The assumptions required for these 
basic equations to hold are too limiting and do not allow 


for all of the processes that are known to occur in combat 
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to be considered in the force attrition model. Some of the 
more common functional forms for attrition rates that have 
been considered and used in Lanchester-combat models are 
shown in Table 2. 

Of special interest is the Tjf functional form. fThis 
functional form is particularly useful in describing attri- 
tion in target rich environments. Peterson {Ref. 8] hypoth- 
esized that the equations, 

dx/fdt = -ax and dy/fdt = -by, 

characterize the early stages of a small unit engagement in 
fich the vulnerability of a force dominates its ability to 
acquire enemy targets. It is recommended that Peterson's 
TIT functional form be considered in the heterogeneous model 
described in Chapter VI to model attrition in the initial 
phase of the battle. The remaining sections of this chapter 
Will discuss several not-so-common functional forms that 
have been formulated and used in "recently" developed combat 


models. 


A. HELMBOLD-TYPE EQUATIONS 
An alternative functional form for homogeneous-force 
attrition rates was proposed by R. Helmbold. dHelmbold 


believed that when opposing force sizes were grossly ine- 


gual, the larger force experienced "inefficiencies of scale" 
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Table 2: FUNCTIONAL FORM NOTATION 


FUNCTIONAL FORM DIFFERENTIAL EQUATION 

Fale dx/fdt = -ay 
dy/fdt = -bx 
PT ler dx/fdt = -axy 
dy/fdt = -bxy 
PI FT dx/fdt = -ay 
dy/fdt = -kxy 
TIT dx/fdt = -ax 
dy/fdt = <-by 

(F+T) | (F+T) dx/fdt = -ay - cx 

dy/fdt = -bx + qy 


in producing casualties. Limitations on space available for 
Maneuver, reduced reaction time, increased command, control, 
and communication problems, etc. result in the reduction of 
casualty producing effectiveness. These limitations nay 
well prevent the larger force from using its full destruc- 
tive capability against the smaller force. There is no 
definition for "grossly" unequal force sizes, but the fol- 
lowing hypothesis is proposed. The validity of using Heln- 
bold's equations in a combat model is dizectly related not 
Only to the force ratios (i.e., xX/y or y/x) but to the rela- 
tive size of the twe cpposing forces with respect to their 


Surrounding terrain. For example, a 5 to 1 force ratio may 


66 





well be a situation where Helmbold's equations apply if both 
of the opposing force sizes are relatively large with res-~ 
pect to the terrain available to them for movement, conceal- 
ment, etc. Conversely, a force with a 30 to 1 force ratio 
May not exhibit inefficiences of scale if the opposing force 
sizes are small in relative sizes and with respect to their 
surrounding terrain. Consequently the decision to use Heln- 
bold's functional form for attrition rates rather than some 
other form (like Lanchester's equation for modern combat) 
really depends on the tactical "Situation". 

Holmbold's functional form introduces a modification 
that accounts for the reduced fire effectiveness of the lar- 
ger force. [In mathematical form his equations would read: 


dxfdt = -a g(x/y) y with x(0) = x, 


dy/dt -b h(y/x) x with y(0) Va 
where a and b represent time-dependent attrition-rate coef- 
ficients, and g(x/y) and h(y/x) represent functions that 
modify fire effectiveness of individual x and y firers 
respectively. According to Helmbold's hypothesis the two 


functions g(x/y) and h(y/x) must satisfy: 


pee = h(t) = 1 


g (u) h (u) E(u) 


E(u) is a strictly increasing function of u. 
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So, the generalized Helmbold-type combat equations can be 


written as: 


dx/dt -a E(x/y) y 


dy/fdt = -b E(y/x) x 
For the case in which E(u) = u° With ¢c a constant 
greater than or egual to zero, the following equations for 


Helmbold-type combat are obtained: 


dx/fdt = ~a ixsy) i y With x(Q) 


it 
”% 


dy/dt = -b (yvxy'™ x with y(0) 


Yo 
Wis commonly called the "Weiss parameter" with W = 1-c. 

The above equations represent the special exponential 
case of Helmbold combat. If a(t) and b(t) represent constant 
attrition rate coefficients denoted by a and b respectively, 
then, when W = 1, the above equations reduce to Lanchester's 
Equations for Modern Combat. When W = 1/2 Lanchester's 


Area~fire Equations are obtained. 


B. HOW TO CHOOSE THE RIGHT FUNCTIONAL FORM 

A very logical question that readers of this study may 
ask at this time is, "Now that I've seen some of the Lan- 
Chester-type attrition rate functional forms, which one 
Should I choose for my analytical model?" The answer to 
this question is that there is no answer to this question. 


The choice of functional form is purely arbitrary and is 
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usually left to the discretion of the model builder. The 
following steps, however, may help in this choice. 
1. Carefully consider the scenario that the combat model 
is intended to portray. The scenario may well indicate 
that some particular subset of the functional forms is 
preferable. 
2. List the assumptions that must be made for each of the 
functional forms in the subset to hold and see which 
assumptions are acceptable to you and your model and which 
are not. 
3. Observe the parameters involved in each functional 
form and eliminate those forms for which the parameter 
values or estimates are unobtainable or for which you are 
not willing to plug in arbitrary values. 
At this point, the final selection of functional form is 
purely one of subjective judgement. 

There is no requirement, however, that a model builder 
must accept one of the already developed functional forms. 
There is no doubt that other forms can be formulated, but 
the formulation will require extensive research, analysis 
and work. A good check to see if an attrition rate equation 
is reasonable is to insures that the Lanchester classical 


equations are special cases of the new functional forn. 
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Once the final functional form has been selected or 
formulated, it may be enhanced as necessary to reflect more 
of the complexity and realism of combat. Some additional 
operational factors that can be considered are; 

Unit breakpoints and battle termination criteria; 

Suppression; 

= ; 

Reinforcements or replacements. 
The main assumption that allows these operational factors to 
be "added" to the attrition rate equation is that synergis- 
tic effects are not considered in Lanchester-type force-on- 
force combat models. No synergistic effects imply that the 
Principle of Superposition holds, therefore, the effects of 


the operational factors are independent and additive. 
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V. FUTURE ENHANCEMENTS 

There are numerous enhancements to the analytical model 
that can be made. The enhancements that are recommended, 
for the most part, are not aimed at improving the inner 
workings of the current model, but at building successive 
models that are more detailed, require fewer limiting 
assumptions, and are generally more complicated. 

For follow-on models that retain the homogeneous force 
concept, it is recommended that the model have the ability 
to allow the user to choose alternative functional forms for 
the force attrition rates. This should be done to better 
describe «he situation of the battle at a particular time. 

A likely method would be to change from one functional form 
to another during the battle in order to more realistically 
portray the attrition occurring then. This enhancement 
would require considerable effort in developing the maxinunm 
likelihood estimates for the various parameters used in each 
alternative functional form, additional work with the STAR 


postprocessor, and Significant change to the current DAMSTAC 


model. Questions which arise from this suggestion are: 


Which form(s) should be chosen? 
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How do you know when to change functional form? 


These questions are partially answered in Section IV-C. 
Further research into this area is required. 

The movement subroutine of the model needs to be modi- 
fied to more realistically portray the movement of units on 
the battlefield. Even though a relatively low-resolution 
time-step operation is occurring, the resolution of the 
movement of units can be refined to vary the rate of move- 
ment from one time step to the next, if necessary. Cur- 
rently, a change in rate of movement occurs only at the 
initial contact between forces and at breakpoints. The rate 
of advance of the RED forces should be allowed to vary based 
on decision logic which considers attrition rate, terrain, 
etc. BLUE forces should be given the capability to move to 
alternate or fall back positions when the battle situation 
dictates. Decision logic must be enhanced to allow BLUE 
forces to move to alternate or fall-back positions. It 
Should be noted that this “pace of battle" question has 
never been fully explained or understood, however, the move- 
ment routine logic in DAMSTAC should parallel that of STAR. 

The concept of units having different formations is 


ignored in the current model. Units are identified by a 
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point location only. Some work needs to be done to allow 
the units to take on different formations for road marches, 
attacks, etc. In order to do this, additional attributes 
such as frontage, depth, orientation and geometric configu- 
ration should be considered to describe the formation of the 
unit. A postprocessor for STAR could take the individual 
tank locations from the STAR data base and determine the 
above mentioned attributes for input into the DAMSTAC model. 
Attrition rates on both sides must be adjusted for the vari- 
ous formations of an attacking unit. 

The inclusion of indirect fire, suppression, smoke, 
intelligence, and close air support are all important to 
making the model a usable analytical tool. Additions of 
communications and logistics, primarily ammunition and 
petroleum, oil, and lubricants (POL), are also important 
enhancements that must be included in the nodel. 

As was seen in Section III-E-3, pa Ne probability that 


a target is visible and acquired in the steady state is 
A n 


Both of the factors represent a particular element of the 


rd 
tt 
rr 
r 
it 


attrition process, one the terrain affects on line of sight, 
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and the other, the acquisition process. It should be appa- 
rent that these two factors can be separated from each other 
with the only common variable, yp , the rate of transfer fron 
the visible state to the invisible state. With this factor- 
ization, if the line of Sight parameters for a particular 
terrain type are known and the attrition rate coefficients 
have been determined, the target acquisition capability of a 
weapon system against another weapon system can be deter- 
Mined. This would allow the target acquisition parameter to 
be used for other terrain types to compute attrition rate 
coefficients without having to estimate the coefficients 
from a high resolution simulation. 

The concept works in this manner. Runs are made on the 
high resolution model (STAR) uSing a particular terrain. 
The attrition coefficients are determined using the maxinuan 
likelihood estimates described in Section III-D. With pre- 
determined data for the terrain coefficients, the line of 
Sight parameter can be factored out of the attrition rate 
coefficient, A. This would leave the value of the acquisi- 
tion factor which should theoretically be constant for the 
firer-target type combination used in the scenario. With 


this value and the terrain coefficient for other types of 


terrain, the attrition rate coefficient can be deternined 
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for a new scenario with the same firer-target combination, 
without running the high resolution model again. This 
method assumes that the terrain affects and the target 
acquisition process are independent of each other. The 
impact of this concept is enormous in terms of time and cost 
saved from not having to make multiple runs of the high 
resolution model. 

One concept which is not strictly an enhancement but is 
an important consideration is the idea of using simulated 
combat data from a source other than STAR. Specifically, 
data generated at the National Training Center from actual 
field exercises could conceivably be used to compute attri- 
tion coefficients in much the same manner as STAR. What 
this provides is not only an alternate source of data, but 
also a different model, since it would include many human 
factors that would not likely be accounted for in STAR. 
Depending on the intended use of the nodel, one source or 
the other may be desirable. 

I* 1S necessary that a progression from homogeneous 
forces to heterogeneous forces be made if DAMSTAC is ever to 
be really useful. This enhancement entails a significant 
amount of research, model building, and statistical analysis 


in order to implement the concept of heterogeneous forces 
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into the analytical model. Target priorities, fire 
allocation algorithms, and a library of attrition rate coef- 
ficients must be developed for various scenarios and force 
mixes. Changing the model to a heterogeneous type requires 
some Changes in the equations (13) and (14), and thus some 


significant modifications to the postprocessor for the 


model. 
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DAMSTAC, in its current configuration, nodels combat 
between two opposing homogeneous forces, but actual combat 
consists of many different weapon-system types operating 
together as combined arms teams. For example, there may be 
infantry, tanks, antitank weapons, artillery, etc. on each 
side. A necessary extension or modification, therefore, to 
DAMSTAC is the consideration of combat between heterogeneous 
forces. This chapter briefly indicates how the basic ideas 
and methodologies used previously to model combat attrition 
for homogeneous forces can be extended and adapted to the 
heterogeneous force structure. 

A key assumption that must be made in order to model 
combat between heterogeneous forces, where Lanchester-type 
differential equations are used, is that the attrition 
effects of various different enemy weapon systems types 
against a particular friendly target type are additive. In 
other words, synergistic effects between weapon system types 
are not considered (Ref. 9jJ. Also inherent in the modelling 
of heterogeneous force combat is the need to introduce an 
allocation factor. This allocation factor may be defined as 


the fraction of the firers of a particular type that engage 
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a particular target type. The allocation factor concept 
will be discussed in greater detail later in this chapter. 
One theoretical concept for modelling the attrition pro- 
cess between two opposing forces is to consider attrition as 
a nonstationary Markov process where the states of the pro- 
cess are defined as the numbers and types of surviving con- 
batants and the ranges between opposing combatants. The 
transition probabilities between the Markov states, under 
this concept, depend on the ratio at which surviving combat- 
ants or weapon systems inflict casualties on a particular 
type of opposing combatant. Thus the nonstationary Markov 
process model provides a means of relating the sequence of 
states and the times of the state transitions observed dur- 
magea battle to attrition rate coefficients. In turn, this 
allows for the estimation of attrition rate coefficients on 
the basis of the observed sequence of states and transition 
times. Andrighetti [ Ref. 10] developed such a nenstationary 
Markov model that related weapon system effectiveness to the 
time sequence of casualties observed in a two-sided, hetero- 
geneous force land combat simulation. There is a flaw, how- 
ever, in Andrighetti's work that eee Should note. de 
implies that his methodology is easily extended to time-de-~ 


pendent at*+rition rate coefficients. This is not true. 
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Andrighetti's thesis does contain a good discussion on 
the theoretical background which allows for the maxinun 
likelihood estimation of parameters used in attrition rate 
equations. An important observation made by Andrighetti ina 
his thesis is that the difference-differential equations 
obtained from the Chapman-Kolmogorov equations of the 
nonstationary Markov model are susceptible to a general 
solution only in special cases as Gordon Clark showed in his 
doctoral dissertation. Clark used the Chapman-Kolmogorov 
equations to develop expressions for time state probabili- 
ties and the expected survivors in a heterogeneous force 
battle. 

The extended version of the DAMSTAC model which shall be 
called DAMSTAC II should be a heterogeneous force analytical 
model that uses STAR output to compute parameter values. 
These parameter values should be allowed to change as the 
distance or range between opposing forces vary during the 
battle. In DAMSTAC, an attempt was made to allow the paraa- 
eter values to change with respect to time. However, it was 
discovered that the maximum liklihood estimators developed 
for the critical parameters in the attrition rate equations 
could not be used for time dependent variables. Two feasi- 


ble alternatives that allowed for variable parameter values 
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were considered. Both alternatives are recommended for 
future evaluation. The first alternative looks at dividing 
the battle into phases based on battle intensity as measured 
by the number of shots divided by the number of live firers. 
(Refer to Section III-E) The second alternative for divid- 
ing the battle into phases is based on the range between 
combatants. Here also the parameters can change from phase 
+o phase. This technique will work assuming that the con- 
batants on each side move aS a Single force and that units 
are located by their centers of mass. For example, a battle 
can be divided into three distinct time intervals or phases 
of combat: an artillery prep phase during which all combat 
vehicles and weapon systems, other than artillery, are not 
participating; a moderate to long range battle phase; and, a 
terminal shorter range battle phase. Parameter values can 
then be computed, using procedures to be discussed later, 
for each battle phase. This second alternative is used in 


the COMANEW model and is the one recommended for DAMSTAC II. 


A. THE ATTRITION-RATE FUNCTIONAL FORM 

The attrition rate functional form that will ultimately 
be chosen for the heterogeneous force model is, of course, 
the prerogative of the modeller. The authors of this the- 


Sis, however, would like to recommend two related attrition 
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rate functional forms for consideration. The first, which 
shall be labeled as Functional Form 1 (FF1), is based on the 
assumption that during the course of a battle fire is 
directed at the highest priority target surviving and 
acquired. Also it is assumed that the relative priority of 
*arget weapon types, as considered by each firer, remains 
constant throughout the engagement and that this relative 
priority is identical for each firing weapon type. Func- 
tional Form 1, developed by Gordon Clark for use in the 


COMAN model, can be written as follows: 


e* (mn) = 0 aie ny pM (2 - P™) 
bj m n) ms ts Aik ny P; Geiecp 
(16) 
ay b a 
tek (mn) % 2 Bikj m, 4; =a = q %) 
where 
a (m,n) = conditional casualty rate of BLUE weapons of 


ype j during phase i of the battle given strengths of a 


and n. 


Cy) 
f+ (m,n) = analog of ty; for RED weapons of type k. 


8 1 





n, = the number of surviving RED weapons of type k 


the number of surviving BLUE weapons of type j 


=] 
it 


Py; = kill rate of a BLUE firer of type j for an acquired 


RED target of type k in the i'fth time interval. 


iy = kill rate for a red-firer/blue-target combination 


corresponding *0 Bi, 


q,;= the probability that a specific RED weapon is unac- 
quired by an individual BLUE firer in the ifth phase of 


the battle. 


pi; = the probability for a specified blue-target/red-firer 


combination corresponding to q;. 


f 


m = <he number of surviving BLUE weapons of higher prior- 


ity than a type j weapon. 
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oe 


n, = the number of surviving RED weapons of higher prior- 
ity than a type k weapon. 

The reader should note that since the assumption is made 
that a firer will engage the highest priority target he can 
acquire, the probability of a BLUE firer engaging a RED 
weapon of type k is g™*(1-q" . Similarly, the probability of 
a RED firer engaging a BLUE weapon of type j is BM (1-p%) . 
These two probabilities may be regarded as the allocation 
maeecors for Functional Form 1. 

All of the critical parameters necessary to compute the 
conditional casualty rates, fy, (m,n) and £". (ten), can be 
obtained either directly or indirectly (using maximum like- 
lihood procedures) from the output of a high resolution 
Simulation such as STAR. 

If one is not willing to assume that target priorities 
remain constant throughout the battle then FF1 will not be 
sufficient. The terms in Clark's functional form which 
result from the assumption of constant target priorities are 
q™ and pi where q®™ is the probability that all surviving RED 
weapons of higher priority than type k have not been 
acquired and pm is the corresponding probability for type j 


blue weapons. 
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One way e extend Clark's formulation in order that the 
fixed priority may be dynamically overridden is by replacing 
the terms qk with 63, where 6}, is defined as the probability 
that BLUE weapon type j engages RED weapon type k, given 
that BLUE weapon type j acquires at least one RED weapon 


type k. Note that 


@(1-g"™ ) = P(j engages kjat least one k acquisition) 


P(acquire at least one k) = P(j engages k) 


Functional Form 2 (FF2) may be written as; 


fy, (ten) = y Oye TM dy (1B; ) (17) 
(at b ny 
fey (m,n) = Do Bin, Dy Oy (1-G;") (18) 


jst 


where Oy; is the RED firer BLUE target probability corre- 
sponding to O, and all other parameters are as defined for 
FF1. Just aS p,q,a,8 are assumed constant over the entire 
phase of the battle, 6 and 6 , for each j and k, are also 


assumed constant over the entire interval. 
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Serious consideration should be given to incorporating 
FF2 or some Similar functional form into DAMSTAC II because 
STAR does allow target priorities to change during a battle. 
Target priorities in STAR are based on the types of targets 
on the firer's detected list and the type of ammunition the 
firer has available for use. 

It may also be useful to use a different functional forn 
during the initial phase of combat where attrition rates are 
more a function of the exposure or vulnerability of targets 
than they ar2 a function of the firer's capability to 
acquire targets. Clark recommended the use of Peterson's 
Logarithmic Law to model attrition during the initial phase 


of the battle. [Ref. 11] 


Be. MAXIMUM LIKELIHOOD ESTIMATORS 

The ability of the DAMSTAC II model to provide insight 
into the interactions represented in the STAR simulation is 
based on the estimation of conditional casualty rates (a 
and B) and other parameters from STAR data. These estinated 
parameters should be constant within a phase of the battle, 
but independent and unrelated to results in neighboring time 
intervals. The independence of parameter values between 
time intervals implies that the estimators for a specific 


time interval, or phase of the battle, can be defined by 
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analyzing a sample of data from the high resolution simula- 
tion observed during that single-time interval. 

The maximum likelihood estimators (MLE's) for the 
unknown parameters a,8, p and g used in FF1 were developed 
by Gordon Clark. It may be of interest to the reader to 


not? sha* simplified versions of these MLE's were used iz 


th 


(0 


DAMSTAC model for homogeneous forces, however, since we 
are now concerned with heterogeneous forces, the "una- 
bridged" versions are presented. The maxinum lixelinood 


estimate for q@ aS a function of p is: 
ol 


= e 
KA 4 ; . 
Meys m 
vD AS (l - Dp Ne Cty = Fy-1 ) 


ato 


} (19) 
9 otherwise 
where a, (pl= the estimator for the parameter q@ for a given p. 


The MLE for p is found by solving: 


Yat 





Mee(D) d. (L - Cy) Cmyay - (fg, + Myg,) 0 Y8Y 
q*) r 
og o(l - 5° tt ) 
(20) 
Vel © A 7: = 
SH by, (pny, oa (Ty - Ty ) ° 
yet jot gat 


m 
Gren -(ngy tm )p 7 )) = 0 


The equations for By; (q) and q are similar. 

The eStimators obtained from the above equations result 
crom che data from a singl2 rteolication of «he combat sinu- 
lation. Mor2 precise estimates are obtainable by using a 


larger sample consisting of a number of observations during 
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the interval in guestion from a number of independent 
replications of the STAR simulation. This modification 
entails an additional summation ()° , where x = the total 
x 
number of replications) in each of the estimating equations. 
The data needed for the computations of the MLE's include; 
y = Pee atl number of casualties 
+t. = elapsed time since the beginning of the interval 
until the y'th casualty 
ty = the end of the battle phase or the elapsed time since 
the beginning of the phase until the battle ended if it is 
less. 
my; = the number of BLUE survivors of type j just prior to 
*he y'th casualty 
Ny, = the number of RED survivors of type k in the battle 
gMSt prior to the y'th casualty 
Cy = 1 i£ the y'th casualty is BLUE ; 0 otherwise 


irer weapon type for the y'th casualty 


rh 
< 
ll 


Jy = target weapon type for the y*'th casualty 


it 


my; = the number of BLUE survivors of higher priority than 
type j just prior to the y'th casualty 

; b 

My; = »s Ayu 


ust 


Ry, = RED analog of my} 
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ahs the number of BLUE weapon-type j casualties due to 
RED weapons of type k 


One technique of computing the MLE*s for @ and D is to 
numerically determine D uSing equation (20) and then substi- 
tute this value into equation (19) for j = 1,2,...,b and for 


k 1,2,..-,t (b = the total number of BLUE weapon types and 


E the total number of RED weapon types). A recommended 
technique for numerically solving for D, the estimate for p, 
can be found in Clark's work on the Dyntacs combat model 
{Ref. 12]. This scenaieae is similar to the one developed 
for the homogensous case discussed in Chapter III. Basi- 
cally the procedure is as follows. Since p represents a 
probability, it must satisfy the inequality O<sp<1. There- 
fore, the function # Ln L({p) can be evaluated at points in 
the interval O<ps1 until two points ,p, and p, , are found 
such that the product of 2. Ln L(p,) and na Ln L(p,) is less 
than or equal to zero. Since 2 Ln L(p) is a continuous 
function, & Ln L(p) = 0 for some p satisfying p, <p < p,. 
This interval is reduced by evaluating 2 Ln L(p) at the 
midpoint of the interval (p,,p,). Repeating this process 
converges to a solution for 2; Ln L(p) and this solution is 


a local maximum. This technique is similar to the bisection 


method for solving nonlinear problems. 
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With respect to Functional Form 2, the introduction of 
the Oj, allocation terms adds considerable complexity to the 
maximization of the likelihood function. One way to circun- 
vent this problem is to compute a "best guess" for the ei, ‘Ss 
from STAR samples (for the given interval). This best guess 
can thenbe substituted into the @, ‘s prior to maximization 
of the likelihood function, thereby reducing the O's to 
constants. The maximization can then be performed, as in 


Clark's original equations, with respect to p, q,a,8 . 


ie CHE "MODEL 

The computational procedures that will constitute 
DAMSTAC II must be developed by the modelier who chooses to 
carry on the work initiated in this thesis. As a starting 
point, however, the techniques used by Clark in his COMAN 
model would be excellent choices. The adequacy of Clark's 
procedures in DAMSTAC II given a particular Scenario and 
force mix, +o predict attrition in a STAR Simulation with 
the same scenario and force mix is questionable and would 
have to be tested. Detailed comparisons and analysis of the 
Output from beth models will be reguired to answer the ques- 
tion of model comparability and modifications to the con- 
puter code and/or computational procedures in DAMSTAC II 


must be accomplished as necessary. The computational cycle 
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employed in COMAN is based on the assumption that the dis- 
tribution times between casualties in the high resolution 
simulation is exponentially distributed. It uses Monte 
Carlo procedures to generate the time until the next casu- 
alty and also to find the casualty weapon type. 

The following steps constitute the computational cycle 
recommended for DAMSTAC II. 

Step 1. The first step in the cycle sets the time equal 
to zero and sets the battle phase, designated by the letter 
meegual to 1. 

Step 2. Next the total conditional casualty rate, 

Ries n) is computed from the equation 

\'(m,n) >> fy,(mn) + oe (mn) for = Uae 2 (21) 

where fy (m,n) and £0) (m,n) are the conditional casualty 
rates in phase i of the battle, given force strength m and 
n, for BLUE weapon type j and RED weapon type k, respec- 
tively. 2! (m,n) is the parameter of the exponential distri- 
bution describing the times to the next casualty. 

Step 3. Now, uSing the exponential distribution with 
parameter xv (m,n), one Monte Carlo's for the time interval 
to the next casualty. Designate this time interval as t'. 
If the next casualty occurs in another time interval, ie. if 


t+ t' > ti where t,;represents the end of the ifth battle 
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phase, then i is increased by 1, the time is set equal to t 
and the program is directed back to step 2 where a'*(m,n) is 
selected as the distribution parameter. If, however, t + t' 


< t; then the next step is to Monte Carlo for the casualty 


~ 


weapon type. 
Step 4. The weapon type for the next casualty is repre- 


sented by the random variables C, and C, which are defined 


as 
C, = j if the casualty weapon is BLUE weapon type j 

O if the casualty is RED 
C. = k if the casualty weapon is RED weapon type k 


O if the casualty is BLUE 
The Monte Carlo procedure uses the conditional casualty 
weapon probability distributions for Cy and C, . These 


probabilities are computed by 


me (0,2) /d(m,n) 
a) 


fr, (M,2)/ Xmn) 


p (Ge jic = 1,m,n) 


p! (C, = KIC 1,mM,n) 


for i= 1,2,..., I and when variable C is set equal to 1 if 
a casualty is determined to occur in interval i from the 
previous steps. The probability that the casualty is of a 
particular type is computed for each BED and BLUE type 


weapon. The sum of these probabilities should equal 1. 
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Step 5. The next step is to Monte Carlo using the Uni- 
form (0,1) distribution to determine the casualty type. An 
example may help in understanding the above procedure. Sup- 
pose there are two types of BLUE weapons and two types of 
RED weapons. Assume the conditional casualty probabilities 
for each weapon are computed using the equations shown pre- 
viously and they are 0.25 for each weapon type, bi, bz, rl, 
r2. If one orders the weapon types in some arbitrary fash- 
ion such as bi, rt, b2, r2, and then accumulates the prob- 
abilities on the interval (0,1), one can asSign an interval 
+o each weapon type. The example would yield the following 


int2rvals: 


mi — (0,0.25) ; 

mee (0. 26,,0.50) ; 
mm —- (0.51,0.75); 
r2 - (0.76,1.00). 


The uniform random number can then specify the casualty 
weapon type. If the random number is 0.60 then the casualty 
weapon type in our example is b2. 

Step 6. After the casualty weapon type has been found 
the next step is to decrement the number of that type weapon 
by 1. Now increase the value of t by t*, the time to the 


next casualty. If * is greater than or equal to ty which 
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represents the ending time of the last phase cf the battle, 
or if either side is annihilated , the battle has ended and 
the program terminates. If t is less than t, then go to 
step 2. 

Observe that the above computational procedures are sto- 
chastic in nature and different replications of an engage-~ 
ment will likely yield similar but different results. If 
this method is chosen for DAMSTAC II, variance reducing 


replications for each engagement should be performed. 


D. THE STAR POSTPROCESSOR 

It is very important that anyone who attempts to develop 
the DAMSTAC II model understand and be able to manipulate 
the STAR simulation model. In the initial attempts to 
extract data from STAR, certain subroutines must be turned 
off. These routines include the Smoke, Engineer, Close~air 
Support, etc. whose effects on parameter estimates used in 
DAMSTAC II have not been thoroughly analyzed. The number of 
types of combatants played in the simulation will have to be 
regulated so that they do not overpower the capabilities of 
the analytical model. It is recommended that no more than 
three different types of weapon systems for each opposing 
force be attempted when initially formulating the heteroge- 


neous force nodel. 


93 





It is important that the finished version of the 
analytical model be as simple and transparent as possible 
and not encumbered with parameter estimating algorithms and 
computations. This consideration plus the amount of data 
that must be extracted from the STAR output require the 
development of a STAR postprocessor program. As an illus- 
tration, if one considers a Simple battle divided into three 
battle phases with only three types of weapon systems on 
each side, the DAMSTAC II model, as it has been described, 
requires the estimation of 114 parameter values and the conm- 
putation of six attrition rates. 

The postprocessor at Appendix E, developed for the homo- 
geneous model, should provide a good foundation for the more 


complicated postprocessor needed for DAMSTAC II. 


94 





The concept of using differential equations to represent 
combat began with Lanchester's classical equations and has 
developed into a very dynamic field of study. The thseoreti- 
cal background of the fitted parameter method has been well 
documented by Clark and Taylor. However, the practical 
application of the theory lacks this same degree of discus- 
sion and documentation. This paper has focussed on some of 
the techniques of applying Clark's theory to an actual high 
resolution, Monte Carlo simulation. The problems are real 
and, as has been seen, the solutions presented are not 
doctrine. It is hoped that further investigation, elabora- 
tion and documentation of this modelling technique will 
occur in the near future. 

The choice of which functional form to use in a particu- 
lar analytical combat model is entirely up to the discretion 
of the model builder. In choosing an attrition rate func- 
tional form *he model builder should consider the scenario 
porzrayed in the model and the assumptions that are compati- 
Ble with it. 

Finally, it is recommended that the entire program for 


*he model be converted to the UCSD Pascal language on the 
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Apple II. The basic structure of Pascal is more conducive 
to combat modelling simulation than the structure of the 
BASIC language. The concept of entities and attributes 
found in SIMSCRIPT II.5 can be closely emulated in Pascal 
while this is not as eaSily accomplished or clearly defined 
in BASIC. 

The modular structure of procedures and functions in 
Pascal is also an advantage in the enhancement process of 
model building. Modellers can use one procedure to repre- 
sent a particular action or phenomenon occurring on the bat- 
tlefield. Any action which is unable to be modelled can be 
left in the program as a dummy procedure for later work. 

Additionally, the Pascal language allows more descrip- 
tive variables to be used than the two-character variables 
used in Applesoft. This would make the program itself more 
transparent and easier to modify and update. 

The DAMSTAC model is clearly in its infancy stage. 


Through continued enhancement, this model can contribute to 


the ongoing work in the field of combat analysis. 
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A = casualty rate for 
AN = angle in radians 


east made by the 


APPENDIX A 


LIST OF VARIABLES 


location and the next coordinating point. 


B = casualty rate for 


BE = the breakpoint for the BLUE force 


BM 


unit strength 


BF (1) 


BP (TI) 


By (I) 


BT (I,J) 


BX (I,J) 


BY (I,J) 


BO (I) 


BO (0) 


B1 = number of coordinating points for BLUE 


line between the current unit 


RED attritting acquired BLUE targets 


measured counter-clockwise from due 


BLUE attritting acquired RED targets 


movement criteria for BLUE units as a fraction of 


number of RED targets being fired at by BLUE unit I 


most recently passed coordinating point for BLUE 


“nae I 

rate of movement for BLUE unit I in M/SEC 

= attrition(casualties) caused by RED unit 
BLUE unit I in one on one battle 

= x-coord of coordinating point J for BLUE 

= y-coord of coordinating point J for BLUE 

current total BLUE force strength for time 


initial BLUE force strength 


oF 


J on 


anit | 
une £5 


period I 


(anetidiac 





imieial Location) 


B1 (I) initial strength of BLUE unit I 
DB(I) = number of casualties or amount of attrition for 
BLUE unit I 

DH = distance travelled by a unit in time increment DT, at 
speed RT 

DI = distance between current unit location and next turn 
point 

DR(J) = number of casualties or amount of attrition for RED 

unit J 

DT = time increment 

DX = x component of distance between current unit location 
and next coordinating point 

DY = y component of distance between current unit location 
and next coordinating point 

D1 = x component of distance moved in time increment DT, at 
speed RR(I) or BR(TI) 

D2 = y component of distance moved in time increment DT, at 
speed RR(¢I) or BR(I) 

M = number of BLUE units 

MT = max time of simulaticn 

N = number of RED units 


N$ = file name (used recursively) 
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P = probability that a BLUE target is not acquired by 
a RED unit 


3.141592654 


PI 


PR = status of printer; 0 printer is off 


1 


printer is on 
Q = probability that a RED target is not acquired by 
a BLUE unit 
RE = the breakpoint for the RED force 
RF(J) = number of BLUE targets being fired at by RED unit J 


range between a BLUE unit and a RED unit 


RG 
RM = movement criteria for RED units as a fraction of 


unit strength 


RP({I) = most recently passed coordinating point for RED 
moiec - 
RR(Iy = rate of movement for RED unit I in M/SEC 
RT (J, I) = attrition (casualties) caused by BLUE unit I on 
RED unit J in one on one battle 
Maqt,J) = x-coord of coordinating point J for RED unit I 


feat, J) = y-coord of coordinating point J for RED unit I 


RO({I) = current total RED force strength for time period I 


RO (9) initial total RED force strength 
R17 = number of coordinating points for RED (including 


initial location) 
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R1(1) 
58 (1) 


SR (T) 


initial strength of 
size of BLUE unit I 


size of RED unit I 


T = simulation time 


XB (I) 
XR (T) 
YB (I) 


me (1) 


current x-coord for 
current x-coord for 
current y-coord for 


GUesen«t Y-GOOrd Lor 


RED unit I 


BLUE unit I 


RED unit I 


BLUE unit I 


RED unit I 
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APPENDIX B 


NOTES 


iz 
\H3 
is 
IO 
to 
rx 
rg 
CS 
ty 
mo 
toe 
1O 
oO 
tem 





A= 2.4 
B= 0.03 
BE = 0.3 
Br = 0 
DT = 30 
MT = 1500 
P= 0.99 


PI = 3. 141592654 


Q= 0.93 


2. Once contact is made (range S$ 2700), the rate of advance 
decreases to 2 meters/second. 

3. The option exists within the program to allow the nove- 
ment of units based on a criteria of fractional force 
strength. When unit strength goes to 50% (movement crite- 
ria), the unit will move to an alternate location (for BLUE 


units) or retreat along its path of advance (for RED units). 
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This appendix provides a brief description of the util- 
i*+y programs associated with the Aggregated Model of STAR on 
+he Apple II Computer. These programs provide the capabil- 
ity to the user to transfer the necessary data from the STAR 
model to this model with the minimum effort. 

The utility programs available are: 

Force Maker 

Results Reader 
Both programs are interactive with the user, who answers 
questions provided by the computer. These programs use the 
commands for the Apple II Disk Operating System (DOS) which 
allow them +o read/write Sequential Text Files and Random 
Access Files onto the disk. Familiarity with these types of 
files and the DOS would improve the understanding of these 
utility programs. [Ref. 13] 

In the following program descriptions, the variable name 
used in the program is listed in parentheses after its 


description. 


102 





A. FORCE MAKER PROGRAM 

This program will write data files containing informa- 
tion about the two oppcsing forces. 

Force Maker produces a Sequential Text Pile containing: 

Force size (number of units) (M) 

The strength of each unit (SB(I)) 

The number of coordinating points for movement of each 
unit (the first coordinating point is always the initial 
location of the unit) (B1(TI)) 

The grid coordinates of each coordinating point for each 
unit (BX(I,J), BY (I,J) ) 

The user has the option to make a file for the BLUE 
force or the RED force. After all the data has been 
entered, a data check routine will print the information for 
all units, so that a check for errors can be made. After 
the user has checked the data and identified the name for 


en 


18) 


file, the program will record the file on the diskette. 
A copy of the program is included as Listing B in Appen- 
dix F and a sample session using this program is attached as 


meugure 12. 
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Be. RESULTS READER PROGRAM 


This program reads the data file of time intervals and 
corresponding BLUE and RED force levels produced by the 
model. Results Reader produces two sequential text files: 
moe first file with suffix ".COEFS" contains: 

Conditional Casualty rate, q@ (A) 

Conditional Casualty rate, B (B) 

Probability that a BLUE target is not acquired by a RED 

firer (P) 

Probability that a RED target is not acquired by a BLUE 

firer (Q) 
The second file with suffix ". RESULTS" contains: 

The number of data points for BLUE force level versus time 

(N) 

The number of data points for RED force level versus time 

(N) 

Time of the battle at time interval I (T) 

BLUE force level corresponding to time interval I (BO(I)) 

Tim2 of the battle at time interval I (T) 

RED force level corresponding to time interval I (R0O(I)) 
mie" .RESULTS" file contains some duplication of data for 
the following reasons. With the ".RESULTS" file written in 


this format, and free of additional data such as that in the 
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" COEFS" file, it can be used directly by the Apple Plot 


Program *o oreduce a plot of tne force levels versus time. 

A listing of «he program is included as Listing C in Appen- 
dix F. A sample plot cf the DAMSTAC output data produced by 
“he Aople Plot Program is included as Figure 13. A sample 
of «h2 outpu* of DAMSTAC from Results Reader is included as 
Figur2 14, 
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Figure 14: 
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THE TERRAIN MODEL 

In addition to the basic computer model, work was done 
+o input the STAR terrain data base into the Apple II compu- 
ter. This effort was done during the initial phase of the 
thesis work based on an erroneous asSumption that the 
authors made concerning the functional form of the attrition 
rates. This assumption was that the terrain of the battle- 
field and line of sight computations must be explicitly 
played in the model. In the given functional forms for the 
attrition rates, line of sight and terrain are accounted for 
by the probability values of p and q. 

Before this assumption was found to be false, four util- 
ity programs, and a model that computed the elevations of 
the units were developed and are included in this document 
for possible future use. If further work is done with a 
model on the Apple II computer, the ability to input the 
STAR terrain into the Apple II may become an important part 
of the model. Examples of models that might use the STAR 
terrain in the Apple II are smail aggregated models with 


functional forms that explicitly play iine of sight, or 


108 





small, high-resolution models that are suitable on the Apple 


II computer. 


A. THE TERRAIN UTILITY PROGRAMS 


The utility programs that allow the STAR terrain to be 


input in the Apple II computer are listed below. 


Map Maker 
Hill Maker 
Woods Maker 
For2st Maker 
These four programs have two portions each: a file maker 
and a file reader. In general, the file maker portion 
allows the user to input a new file or update/change an old 
file. The file reader allows the user to print out on the 
screen or to a printer the selected file. 
Por the Apple II combat model, a STAR terrain file (eg. 
FULTERR file) is divided into its four main parts: 
the list of hill reference numbers in each grid square; 
the hill parameters; 
the list of forest reference numbers in each grid square; 
the forest ellipse fitting parameters. 
Each part is created, respectively, by one of the four util- 


ity programs: Map Maker, Hill Maker, Woods Maker, and Forest 
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Maker. When these parts are being created on the Apple 
Disk, the same general file name must be used to identify 
the data's terrain (eg. FULDA). Each of the utility pro- 
grams will add its own suffix to the general file name. For 
example, for the Fulda terrain, the four utility programs 
will produce files named: 

FULDA. MAP 

FULDA.HILLS 

FULDA. WOODS 

FULDA.FORESTS 
This suffix identifies the file as a particular portion of 
the terrain data, and will allow the Apple computer to 
recognize it from the other files with the same terrain name 
when reading the file for the combat model. 

1. Map Maker Program 
This program will read and write data files contain- 

ing a catalog by grid square of the hills which influence 
the *errain in that grid square. This program produces an 
equivalent file to that in STAR called LIST.H. Map Maker 
produces a Sequential Text File containing: 

The size of the battlefield (normally 10 by 10) (L,W) 

The lower left grid square in 1000 meter units (LX,LY) 


The base elevation of the battlefield (HL (0,0,0)) 
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For each grid square I,J: 
(1) The number of hills influencing the terrain 
(HL (I,J, 9) ) 
(2) The reference number of these hills (HL(I,JdJ,)) 
The user has the option to make a new file or access 
an old file for updating or changes. For new files, the 
program will start with the lower ieft grid square and pro- 
ceed from west to east then south to north. The data can be 
saved after the information for a grid square has been 
entered. For old files, the program will access it, print 
out the file on the screen and ask for the starting grid 
Square. The grid square index is used to identify in which 
grid square to start. The grid square index 1S an integer 
from 1 to Lin the east-west direction, and 1 to Win the 
south-north direction of an L by W sized battlefield. The 
lower left grid square is 1,1; the lower right is L,1; the 
upper left is 1,W; and the upper right is L,W. 
A copy of the program is attached as Listing D in 


Appendix F and a sample file is attached as Figure 15. 
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This program reads/writes data files containing the 
fitting parameters for each hill on the battlefield. Hill 
Maker produces a Random Access Text File containing: 

The number of hills on the battlefield (N) 
mor @ach hill Tf: 

(1) X coord of the center location of the hill (H(I,1)) 

(2) Y coord of the center Location of the hill (H(I,1)) 

(3) The elevation of the hill top measured from 
zero=sea level (H(I,3)) 

(4) The orientation angle of the ellipse measured in 
degrees counter-clockwise from east to the major axis. 
(H(I,4)) 

(5) The eccentricity, defined as the ratio of major 
axis length to minor axis length (H(I,5)) 

(6) The spread, defined as the distance in meters meas- 
ured along the major axis from hill center to the contour 
line which is 50 meters down from the peak (H(I,5)) 

(7) The maximum height of the normal curve describing 
this hill mass (H(I,7)) 

(8) The vertical distance measured from the peak beyond 


which this hili is not considered in computations (H(I,8)) 





A copy of the program is included as Listing E in 
Appendix FP, and a sample file is attached as Figure 16. 
3. Woods Maker Program 
This program will read and write data files contain- 
ing a catalog by grid square of the woods which influence 
the terrain in that grid square. This program produces an 
equivalent file to that of Map Maker. Woods Maker produces 
a Sequential Text File containing: 
The size of the battlefield (normally 10 by 10) (L,W) 
The lower lef*+ grid sguare in 1000 meter units (LX, LY) 
For each grid square I,J: 
(1) The number of forests influencing the terrain 
(FL (I,J,9)) 
(2) The reference number of these forests (FL(I,J,i)) 
Woods Maker contains the same options as those found 
in Map Maker. A copy of the program is attached as Listing 
F in Appendix F, and a sample file is attached as Figure 17. 
4, Forest Maker Program 
This program reads/writes data files containing the 
fitting parameters for each forest on the battlefield. For- 
est Maker produces a Random Access Text File containing: 
The number of forests on the battlefield (N) 


For each forest I: 





eRUPREO 
FULDR HILL BATA 


NUMBER OF HILLS: 159 

‘ S32 350 425 142 2 {ez 252 
2 540 953 432 114 z 320 2s 
z 540 352 420 45 5 500 35 
és S42 925 420 74 2 320 is 
5 Shu 954 425 iS x 52a 259 
= S43 345 452 2 i [20 70a 
7 559 345 450 2 1 700 252 
2 Sez 339 450 255 3 322 S20 
3 Sas 1925 471 43 2 =2a 42a 
12 [25 1314 330 48 ; 7200 142 
rt S75 362 462 42 3 3200 210 
iz S44 1Oz2 3S 25 2 422 700 
ie 532 1 O27 325 52 1500 300 
a 572 i213 415 eS z 1220 Ty 
is S32 1025 330 Lé z 320 25 
16 535 1Diz 475 jure 2 300 225 
7 562 1209 330 a e 1720 720 
13 524 1202: 462 {e7 { 522 235 
13 S73 1002 335 12 A 52a 150 
22 So4 237 43S 25 : 35 15 
ag S32 S4 430 { { SZ 249 
ae 529 320 z55 16S. A 2200 142 
oy =Si 1220 394 142 Ee 1750 1a 
Zé =34 1205 394 145 é 120g 120 
25 S24 343 = 220 4 2900 142 
25 39 342 240 165 3 2020 122 
ay 52a 245 730 a7 es 1220 14a 
“3 S4F 3252 332 =5@ g 120 15g 
zs S45 354 730 33 a 1220 is 
2 =45 229 TE4 2 : 1 22a Lae 
= 55 351 400 210 = i 20a 15a 
a 530 he s3a 125 e sau £92 
a. =30 355 SEQ 35 3 a3 310 
74 536 353 sa5 35 = 1500 700 
35 S25 352 s25 12 = i 5og 7G2 
“3 Zit S72 a, 2 : [og 702 
- S464 372 279 74a - 750 129 
ze 539 375 740 L20 zs 50a 3 
73 STz 375 472 250 a soa 329 
4 =77 375 | 479 745 2 =20 ace 


Figure 16: Sample File for dill Maker Progran 
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(1) X coord of the center location of the forest 
(F(I,1)) 

(2) ¥ coord of the center location of the forest 
(F(I,1)) 

(3) The height of the trees in the forest (F(I,3)) 

(4) The orientation angle of the ellipse describing the 
forest measured in degrees counter-clockwise from east to 
the major axis (F(I,4)) 

(5) The Length of the semi-major axis of the ellipse 
(F(I,5)) 

(6) The length of the semi-minor axis of the ellipse 
(FPF (I,6)) 

Forest Maker contains the same options as those 
found in Hill Maker. A copy of the program is included as 
Listing G in Appendix F, and a sample file is attached as 


Figure 18. 


B. THE TERRAIN MODEL PROGRAM 

A copy of the model which explicitly plays terrain is 
included in Appendix F as Listing H. A line drawn along the 
edge of the program delineates the terrain model portion 
which was added to the basic model. The additional varia- 


bles required to run the terrain model are also included at 
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the end of this appendix. Professor James K. Hartman [Ref. 
14} provides the details of the methodology for line of 
sight and elevation using the STAR terrain. 
equivalent STAR variables are in ( ) 
Al = computing parameter for a hill (A) 
B1 = computing parameter for a hill (B) 
CR = critical value parameter for a hill (CRIT.H) 
C1 = conversion of orientation angle to radians 
C2 = (spread of a hill) squared 
FI = value of the elevation of a point due to a hill (FTI) 
H(I,J) = the J-th parameter for hill I; the parameters are: 
1 = x-coord of the center of the hill (XC.H) 
2 = y-coord of the center of the hill (YC.H) 
3 = peak elevation of the hill (PEAK.H) 
4 = orientation angle of the hill (ANG. H) 
5 = eccentricity of the hill (ECC.H) 
6 = spread of the hill (SPRD.H) 
7 = max height of the hill (HT.H) 
8 = cut height of the hill (CUT.H) 
H(I,90) = status of hill I; 


0 


hill I has not been checked for LOS 


1 hill I has been checked for LOS 
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the number of the K-th hill in grid square I,J 


ul 


HE) (1 , J , K) 
which influences the terrain in that grid 
Square IT = 1 to L; J= 1to W; K = 1 
to HL(I,J,0) 


the total number of hills in grid square I,J 


He {i,J, 0) 
which influence the terrain in that grid square 


HL (0,0,0) 


the base elevation of the terrain map 

H2 = temporary variable identifying the hill being checked 
For line of sight 

L = number of grid squares on the terrain map in the 

x direction (NGRIDX) 

LX = x-coord of lower left grid square of terrain map in 
1000 meter units (X.LO.BDRY) 

LY = y-coord of lower left grid sguare of terrain map in 
1000 meter units (Y.LO.BDRY) 

M$ = name of files containing terrain data 


NH total number of hills on the terrain map (NHILLS) 


P1 = computing parameter for a hill (PXX.H) 
P2 = computing parameter for a hill (PXY.H) 
P3 = computing parameter for a hill (PYY.H) 
QI = quadratic function of the ellipse represnting the 


hill (QT) 
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x 
i 


number of grid squares on the terrain map in the Y 


direction (NGRIDY) 


os 
it 


current x-coord of unit whose elevation is being 


computed (in 100 meter units) (X) 


XB (T) 
XR (I) 


.> = 


at = 


X2 = 


= current x-coord for BLUE unit I 


current x-coord for RED unit I 


distance in x direction between unit location and 
center of hill mass (XS) 

relative number of grid squares in x direction of a 
BLUE unit location with respect to the lower left 
grid square 

relative numberc of grid squares in x direction of a 
RED unit location with respect to the lower left 
grid square 

number of grid squares in x direction between a BLUE 


unit and a RED unit 


X4,Y4 = current grid square being checked for line of 


Sight in LOS subroutine 


Y = current y coord of unit whose elevation is being 


computed (in 100 meter units) (Y) 


YB (TI) 


YR (I) 


current y-cocrd for BLUE unit I 


current y-coord for RED unit I 
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is = 


Y1 


Y2 = 


Y3 = 


distance in y direction between unit location and 
center of hill mass (YS) 

relative number of grid squares in y direction of a 
BLUE unit location with respect to the lower left 
grid square 

relative number of grid squares in y direction of a 
RED unit location with respect to the lower left 
grid square 

number of grid squares in y direction between a BLUE 


unit and a RED unit 


Z = temporary variable for the elevation at a point ona 


mead (2) 


ZB (I) 


ZR (J) 


elevation of BLUE unit I 


elevation of RED unit J 
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The following appendix is a summary of the work done by 
other students as a part of a joint project for OA 4655. 

The authors gratefully acknowledge the work done by Captains 
Ambrose R. Hock and Steven L. Maddox in developing this 
postprocessor to reduce the amount of work required in 
extracting the necessary data from STAR. Their work is 
included here in this appendix since it is an integral part 
of the parameter estimation process. 

The postprocessor developed for this thesis is a Sin- 
script language program that provides to the user a Summary 
of the data generated from the Simulation of Tactical Aiter- 
native Responses(STAR) Model. Additionally, the postproces- 
sor performs some routines that assist the user in the data 
analysis of th2 model. The STAR Model is a Simscript lan- 
guage program that simulates combat between two combined 
arms teams in a combat environment that includes field 
artillery fire and electronic warfare. 

Since the model assumed in the project oniy required 


data generated by a tank vs tank battle, it was necessary to 
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change certain aspects of the STAR model. Thus it was 
necessary to "turn off" the field artillery and Electronic 
Warfare modules within the STAR model. Also the non-tank 
units needed to be deleted. Referring to figure 19, the 
Electonic Warfare code(EWCODE) was deleted through the use 
of a comment card, anda Electonic Warfare off (EWOFF) 
subroutine was inserted. 


//STARTHES JOB (3102,0234) ,*COMB ARMS BATTLE',CLASS=C 
fesea. SYSPRINT DD DUMMY 
//SIM.~SYSIN DD DISP=SHR,DSN=M4SS.S3102. THESIS (WRKPREAM) 


// DD DISP=SHR, DSN=HSS.S3102. THESIS (SNAPDOTR) 
//*® DD DISP=SHR, DSN=MSS.S3102. THESIS (EWCODE) 
// DD DISP=SHR, DSN=MSS.S3102. THESIS (MOVCOORD) 
// DD DISP=SHR, DSN=MSS.S3102. THESIS (COMMWORK) 
// DD DISP=SHR, DSN=MSS.S3102. THESIS (TEMPCOMM) 
// DD DISP=SHR,DSN=MSS.S3102. THESIS (TEMPGRND) 
7/ DD DISP=SHR, DSN=MSS.S3102. THESIS (EWOFF) 

Vf DD DISP=SHR, DSN=MSS.S3102. THESIS (FAEXPER) 
i / DD DISP=SHR, DSN=MSS.S3102.THESIS (ARTYMESG) 


Figure 19: Electronic Warfare Deletion 


Thus, whenever the main program called for the usage of EW, 
the EWOFF subroutine automatically returned the program back 
to the point where the EW routine was called and the sinula- 
tion continued with no Electronic Warfare generated. The 


field artillery within the STAR model was cancelled by 
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having the Simulation create FA fire at time equal to 5000 
time units(see Figure 20). 


Schedule a FABEGIN in 5000.0001 time units 


Pigure 20: Field Artillery Deletion 


Since the Simulation was scheduled to run only 2500 time 
units, no field artillery was generated. Finally, the 
internal logic of the STAR model was used to "destroy" the 
non-tank units within the simulation by changing their basic 
load to zero(see Figure 21). This action created situa- 
tions, where the tank elements did not consider the non-tank 


elements as lethal or dangerous targets and thus were 


ignored. 
7 1 1 190 38 1 #10 17 110 100 00 00 
ze 3 0 00 Q 0 00 0 000 0 00 00 
2 4 1 QO 1 10 0 00 0 000 0 0000 
a 6 0 00 0 10: 1.2 000 0 00 00 
q 4 50 20 Q 50 60 30 QO 50 000 0 00 QQ 
7 inet. 022 1 #10 18 110 100 00 00 
2 38 10 14 0 000 110 11 0000 
1 1 1 10 38 1 Ome 1 10 100 0000 
Z 3 0 00 0 0 9000 090 0 0 00 0 0 
2 4 0 000 0 00 0 000 0 OO = Oy 0 
| 6 0 000 0 QQ 90 000 Q 00 0 0 
4 1 50 20 0 50 60 30 0 50 000 0 00 00 
17 Weer 0° 22 1 1090 18 1 1 0 100 Om07 0 0 
ge 3 0 00 0 0 Q0 0 000 0 0000 


Figure 21: Alteration of Basic Load 


Output from the STAR model can be divided into three 


categories: 
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Update list of attacker/defender status. 


Shot data 


Significant event listing (1.€. movement, EW and FA 

events) 
Only the shot data from the STAR model is needed to provide 
the necessay data for computations. Figure 22 is a sample 
of this output. In order to utilize this data, the internal 
output control of the model was used to place the shot data 
in a mass storage location within the computer to be later 
Manipulated by the postprocessor. 

The flow chart for the postprocessor is contained in 
Figure 23; a copy of the program itself can be found in 
Appendix F as Listing I. The program flow is very basic and 
needs little explanation. From the mass storage location of 
the computer, the postprocessor reads pertinent information 
into the SHOTLIST. Once all of the data has been read, the 
program processes the data into a casualty list(CASLIST) 
keying from a status of "dead" from the Shotlist. Figure 24 
is a sample of the casualty list output. Once CASLIST is 
created the postprocessor performs histogram computations 
and graphs. The postprocessor can print two types of 
histograms. Figure 25 contains the Simcsript Library pro- 


gram for the histogram. Additionally, the HISTG routine in 
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FORTRAN can be accessed to produce a histogram of the data. 
After the histogram data is generated, the postprocessor 
calculates the attrition rate coefficients using the formu- 
las discussed in Section III-E of this thesis. 
of the postprocessor (see this appen- 
dix for variable description): 

Preamble: 

lines 1- 9 definition of permanent entities. 


lines 10 - 39 definition of temporary entities. 


lines 40 = 54 defanition of data to be collected for 


histograms. 
Main: 
line 2 variable definition. 
line 3- 4 system parameters read. 
lines 13 - 31 creation of shotlist, reading data from 
tele. 


lines 43 = 89 creation of casualty list, file "dead", 
lines 90 - 125 processing of histogram data 
lines 125 - 186 processing of attrition rate coeffi- 


cients. 
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DEFINE 
VARIABLES 














READ 
DATA 
No 
CREATE COUNT 
> PROCESS 
SHOTLIST eaten 
COMPUTE 
HISTOGRAM 
INITIALIZE PROVAL 
DATA = 
ak 
FOR P=.OL 
To . SY .O! NO 
INITIAL/ZE 
Dyas 
y 


Figure 23: Flew Chart of the Postprocessor 


OAS, 









at 


> 
NTABLE 


CALCULATE 
¥ AND 2 


FOR EACH 
CAS IN 
BASES T 
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NO 


— Se 


PRINT 
OUT PU T 


SHOT LIST DATA FOR CASUALTIES FOR THIS STAR SIMULATION 


NUM Tiae ae a Sa 2YPE PERE BETWEEN Sarees 


1 7 

2 334 46 1-7 12 

3 468 58 US ay 134 

dh 574 28 Va 106 

5 58 1 42 7 7 

6 591 27 US] 10 
107 1210 72 1-7 6 
108 1238 69 1=7 28 
109 1265 73 aes 259 | 
mo 861303 is ee), 38 
it 1323 Us 1-1 20 
112 «1421 76 i=) 98 


Figure 24; Casualty List 
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NOMBER OF BLUE KILLED = 28 


BLYUU = 47.25000 BLSIG = 299.259 
NUMBER OF RED KILLED = 84 
RDMU = 13.24096 RDSIG = 219.106 


HISTOGRAM OF BLUE AND RED CASULTIES 
(DATA BASED ON TIME BETWEEN BLUE CASULTIES 
AND TIME BETWEEN RED CASULTIES) 


INTERVAL=10 UNITS #BLUE CAS #RED CAS 


1 3 11 
2 4 8 
3 4 20 
4 0 4 
5 0 8 
6 0 4 
7 0 0 
8 0 8 
9 0 0 
10 0 ) 
11 0 uy 
12 4 8 
13 0 4 
14 4 4 
15 0 0 
16 0 0 
17 Q 0 
18 0 0 
19 4 0 
20 0 0 
21 5 0 


Figure 25: Simscript Histogram 


RAHAT = estimate of the attrition coefficient for RED 
attriting BLUE 


estimate of the attrition coefficient for BLUE 


BHAT 
attriting RED 
BLHST = Simscript routine for the histogram of BLUE tine 
between casualties 


BLMU = Simscript routine for calculating the mean time 
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between casualties 


BLSHTHST = Simscript routine for the histogram of BLUE 


shot times 


BLSHTIG = Simscript routine for calculation of the 
standard deviation of BLUE shot times 

BLSHTMU = Simscript routine for calculation of the mean 
of BLUE shot times 

BLSIG = Simscript routine for calculation of the 


standard deviation of BLUE time between casualties 


BLUE.SHOTS = time of the previous BLUE shot 


BLUE.TIME = time of the previous BLUE casualty 


CAS = 


casualty 


CASHST = Simscript routine for the histogram of 


CASIG 


total time between casualties 
= Simscript routine for calculation of the 
Standard deviation of the total time between 


casualties 


CASLIST = array containing list of casualties 


CASMU 


CBLUE 


CK = 


Simscript routine for calcualtion of the mean 
for the total time between casualties 

= number of surviving BLUE force elements 
integer used in the calculation of the attrition 


coefficient 
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= 
tl 


BLUE casualty 


© 
tt 


RED casualty 

CRED = number of surviving RED force elements 
CURTIME = time of previous casualty 

FNAME = firer's name 


FSWTYPE = firer's system-weapon type 


J.BL naumber of BLUE caSualties(used in BLHST routine) 

J-RD = number of RED casualties(used in RDHST routine) 

JBLHST = number of BLUE shots(used in BLSHTHST routine) 

JCAS = number of casualties(used in CASHST routine) 

JRDSHT = number of RED shots(used in RDSHTHST routine) 

JSHOT = number of shots( used in SHOTHST routine) 

JTOT = number of casualties(used in TOTHST routine) 

LIM = array of time periods that represent phases of 
the battle 

LIMA = lower time limit of phase 

LIMB = upper limit of phase 


partial derivative of p 


LP 


LQ partial derivative of q 
MAXBLUE = initial BLUE force size 
MAXRED = initial RED force size 


MINP = minimum of p 


MINQ = minimum of gq 
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MK1 = number of RED elements in time interval 

NK1 = number of BLUE elements in the _time interval 
NLIM = number of time intervals to be evaluated 
NSHOT = number of shots in the total battle 

NTABLE = number of tables(battles) to be evaluated 
NOM = counter fo number of casualties 


NUMBL.CAS = number of BLUE casualties 


NUMRD.CAS = number of RED casualties 
NUMCAS = casualty entity number 

P = probability of not acqiring a target 
PHAT = estimated p 


PROVAL = process variable 


Q stop 


— 
i 


process 


QHAT = estimate of q 


ye 0] 
to» 
= 
Q 
td 
il 


Tange of firer to target 


RDHST = Simscript routine for the histogram of RED 
casualties 
RDMU = Simscript routine for calulation of the average 
time between RED casualties 
RDSHTHST = Simscript routine for the histgram of total 
red shots 


RDSHTMU = Simscript routine for calculation of the average 
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red shot time 
RDSHTSIG = Simscript routine for the caiculation of the 
standard deviation of the RED shot time 
RDSIG = SimScript routine for calculation of the 
Standard deviation for RED time between casualties 
RED.SHOTS = time of previous RED shot 
RED.TIME = time of previous RED casualty 
RUNNO = run number 
SHOT = shot 
SHOTHST = Simscript routine for total shots fired 
SHOTLIST = array of shot times 
SHOTMU = Simscript routine for calculation cf the average 
+ime of shots 
SHOTSIG = Simscript routine for calculation of the 
standard deviation of total shots 
STATUS = status of shot 
STATUS.CAS = status of casualty 


T.BTWN.CAS = time between casualty 


im CAS casualty number 


number of BLUE casualties within the interval 


TBLUE 
TGT.CAS = target casualty 
TIME = time of shot 


TIMECAS = time of casualty 
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TIMELIM = maxtime of battle when a casualty occurred 

TMCAS = time of previous casualty 

TMRDCAS = time of previous RED casualty 

TOTAL.CASTIME = total casualty time 

TOTAL.SHOT.TIME = total shot time 

TOTAL.TIME = total casualty time for the creation of 

T.BIWN.CAS 

TOTHST = Simscript routine for the histrograam of 
total casualty time 

TOTMU = Simscript routine for calculation of the 
average casualty time 

TOTSIG = Simscript routine for calculation of the 
Standard deviation of the total casualty time 


*otal RED casualties within the time interval 


TRED 
TSUB = time counter 
TSWTYPE = target system-weapon type 


TTBLUE = «otal BLUE casualties for XX battles 


TTRED total RED casualties for XX battles 
TTCAS = total casualties for XX bat¢les 
TTT = number of shot 


TYPE.CAS = system weapon type of the casualty 


VAR = attrition coefficient variables 


137 








DAMSTAC MODEL PROGRAM 


14 REN WRPOATED 24 FES Be 

24 FEM OETERMINTSTIC S6isREGSTEN HODEL mF STOR ON THE 
21 REM SPPLE COMPUTER © CRMSTHIC ® 

34 QIK mB 5 a Pant 5 a ere a] 1, SRE a ‘te B VAC 5 B a5) I cys a , tH a PY Sean ed r) 1h 3 
44 OM BPOS) FPC52,58¢59,5R0355 Ye SL 

38 OIM BFC S>,.RF6S°%,.8T65,.59,.RT6 5,559,080 57,.0R¢ 5 3 
ba OM BRS >,RR SS) 

‘OQ G“IM B1652,81055 


34 T = 4:0T = 2Q3 MT = {5Ha 
a 0S = CHRS (45 

1646 P = H.3819 = 4,595 

116 PI = Z.1S1592685464:PR = 4 
124 A = 2.4:8 = 1.2 

134 BE = ALS:RE = 45 

itn EM = .SsFM = LS 

154@ KK = 

164 NK = NT + OT 


iSt FEM ZS=FLAG FOR TYPE OUTPHLIT . 
Ibe REM cé=1 SIVES UNIT 

162 FEN LAOCARTIONS & SICES 

ifs FEM fe=4 SIVES TINE & TOTAL 
i6bS REM FED * BLUE FORCE LEVELS 
bb cl = 4 

1PM DIM BAONK OA NK 3 


lel REM FEAT OATA FROM FILE FOR 
lsc BEM S&LUE FREIE. 

133 SEM ORTA INCLUOES: 

134 FEM N=MUIMGER OP UHITS 

135 FEM SL=NUMBER OF COORO PTS 
ts6h OREM OBK,.BY=GRID OF COORO PT 


3 INPUT “NeME DF FILE =OR SLile O9TAT & INg 
2y4 PRINT 05;"NPEN “sGNss". ALOE" 

zl PRINT O3s"RES80 “GNss", BLUE" 

220 INPUT 4 

3m FOR I 
oy sha 


IB ty 


4 
map & 


it 


SoM FOR J = tf To &1 

2hu alle Bee l.dce INPUT Aye D,.d% 
Sr NEST > NEST 

SM PRIWT O3s°CLOSE “shes. Sue 
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es 
252 
233 
aad 
234 
shy 
314 
424 
me 1 
RvR 
eS 
kt 


ool 


re | y ye . 
g 


we ee See ae Oe Oe Oe ae Oe a 
keh Cd Ra = = & ‘Ys Ly ‘Js wn E Hy fa n+ ho 
De Se Se ls Se Ss oH 


oy LY aie 


oy ey A 
HWM F 
See ey & ei Seah 


~- 
7, 


be Ti 


REM ‘SET [NITIRL UNIT 
RFEH LOCATION EGUeL Th 
PEM 1ST CORO PT. 

REM 


POR T = t TO Mex I> = BKC TLR Do = Bee Tt 2 SPeL 


= 1 


PRINT B1: FOR Jt = 1 TO Biz PRINT BME Loe Tete: NEXT J 


PRINT “SIZE OF BLUE UNIT “sia: TNPLIT SBT 
Biel = Sec: 

FEM 

REM ‘SUM INITIAL FORCE STZE 

REM 

cui = Ber as + SBE T 3 


REM RATE OF BLUE MOUVEMENT 
REM EGLIBLIS 


Wea 


REM REAL QATA FROM FILE FOR 
REM RED FORCE. 

FEM ‘SHME FORMAT 48S BLUE 

FEM 

INPUT “MBME OF FILE FOR RED ORTQ? "sNS 
PRINT O3;5"0RPEN “sNS5". RED" 
PRINT O3;"RES0 “sNS:". RED" 
INFUT N 

FOR I = 1 TO 

INPUT Ri 

FMR J = 1 TO Fl 

INPUT Ree [wos INPUT Ree T,.d 
NEST =: NEXT 

PRINT OS;"CLOSE “GNSs". RED" 


POR IT = 1 Th Ns¥RO TS = RXCT.1 oR Ts = Be Tot ts 


PRINT “SIZE OF RED UNIT “sles [PHT ‘Ske 3 
Ricly See I 3 

PAC PCH + SRE I 2 

Ree I 2 

NEST I 


PRINT =: INPUT “Om ‘vot WENT THE PRINTER CN? © Ye 


IF ws = "Y" THEN PR = { 

HIME : UTHR 3 

IF PR = 1 THEN PR# I 

PRINT 

REM 

PEM FRINT [INITIO OATS 

FEM 

PRINT “TIME“s: POKE 36.7: PRINT "SILI": : 


Sik e 


3K ects 


PRIHT 


it CFE ub 





Si 


6 Oe Os HC 


ws ys ss «fi ca An ef) s_ 5h 
fea we Qe Of 
Samah TiN & 


Se ee? es ee We ope ope Fe Js y Re 
(Ain cH esi 4. a a fo hs ho ho 
SDA SSNS OOS 


dD. 


~] >. 


Jee 
ff iS Si) rv gt 


oye te eo 


PRH +4 
REM 


FEM WHICH OUTPUT TYPE? 


REM 


ie 22 3 


1 THEN 
IN-vSUB 143M 


MISA LPRRe BATT Sie 


REM MHIN PROGRAM 


T=T+0 
FEM 


+ 


REM CHECK [FP MAX TIME 


Rew b> Ee CEeseL 


REM 
IF T < 


S056 VF 


BUSHB 11a 


FEM 


= MT THEN 53a 
PRINT “ENO OF BATTLE DUE TO TIME’: BOSUIR 1544: END 


bys 


REM 


>_ /: « > 


MOVEMENT ‘SHIRRGUITINE 


FEM ATTRITION SUBRMITTINE 


seh GHEREMENT THE | CME 


FEM FERIOO 


KEM 
KK = KK + 
BHC KK) = 
REM 


1 


Me Rit KK oo 


TNOES 


=H 


men Sen CUE EOPEE EVEL 
fem PER THIS TIRE PERIOD 


SEM 

FOR I = 
BatKK = 
MEST I 
FEM 


1 Th ™ 


RickKK o + ‘SBC 0% 


REM ‘SUM RED FORCE LEVEL 
REM FIR THIS TIME PERIOO 


FEM 

FOR J = 
RICKK > = 
MEST o 
REHM 


1 


TH ON 


PRCKK oO + SRC) 


FEM WHICH TYPE CIF MITPHITsS 


Rer 
IF 22 = 
BOG 14 


t 
th 


THEN 
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we) ops 
iS 


es CAO Oe ie os Ot es Ce oes 


fej fi. ft ped pt rt NAN 
Bs oe gop Uh gf Gly Ts 


Gis 


oie OT = eRely = COS ¢ 


ee 

mie HECK TF ENC OF BATTIE 
MEM - CRITERIA FOR PORTE 

EM LEWEL'S TS EXCEEDED 

FEM IF CRITERIA [35 ESCEEDED 
FEM GOTO SBVE FESULTS 

REM 

IF BUC KK> > BE * RuC A> THEN Prr 
PRINT : PRINT “END OF BATTLE QUE TO ATTRITION OF BLUE” 
S03UB 1544: END 

IF PUCKK) > RE * Fatdw> THEN FSR 
PRINT =: PRINT “ENO OF BATTLE QUE TO STTSITICN OF RET" 
IIISUB 1See: END 

HiT 52a 

FEM MOVEMENT SUBROUTINE 


REM 

REM FOR FED UNITS 

REM 

FOR T = 1 TON 

REM 

FEM GET LAST COonRPO PT 
REM 


J = RPCT) 


REM 

FEM COMPUTE COMPONENTS OF 
REM OISTGNCE BETWEEN NEXT 
FEN CUOnRO PT GNA CURRENT 
REM LOCATION 

FEM 


CY = RYCT.I & 12m ROT 
OX = RKCLT,I + 13 -— KRCTS 


REM 
FEM OU ISEOMETRY TO INSET 
REM (CORRECT ANGLE «BN > 


REM 

IF 0% = 4 GNO OY * & THEN GN = PT © 7: GOTO aaH 

IF OX = 4 BND OY < @ THEN GN = 2 <= PT ¢ 3: NTH B6H 
BN = HIN (DY - O42 


TF O8 <= & THEN BN = PI + BM 


REM COMPUTE OH, OISTANCE 
REM TO BE TRAVELLED IN 
pete i > Uline Perris 


SIN ¢ 


+ fist a 


2 = OT + RRCT) 
= SHR COL « 


es comes} 


14 1 





Shad FEN 

35 FEM COMPUTE OI. OISTANCE 
346 FEM TO NEXT WOORO PT 

347 FEM FROM CURRENT LOCATION 
343 REM 

314 OT = SAR COX « 2 + OY ww Sh 
311 REH 

412 FEM IF OF < OM, MOWE UNIT 
S15 FEM ONL'Y THE BMDUNT TY 
414 FEM NEXT COORD PT 

315 REM IF OI [5 USED. URPOATE 
316 FEN LAST COORD PT 

StF REM 

3x20 [F QI > = OH THEN 354 

324 O01 = OK:02 = O'r 

‘agua RPC T > = J + 1 

954 YRoT> = TNT cYRC TT + 2% 
S64 HROID = INT CBRE To & O19 
37H NEXT 

474 REM 

475 FEM OD THE SHAME FOR BLE 
3°7B REHM 

tsa FOR [— = 1 TOM 

385 IF BRcI> = & THEN 1158 

394 J = SPC I: 

lyn O'P = Brel. + 13 — ECT» 


1Miw OK = BKC Te, + 13 = SBC TS 

itze IF 0% = 9 SNO OY = @ THEN ON PY “4 te ANITM thBe 
iwsy TF OX = 4 SNO GY < 8 THEN BN a * PT “ 3: GOTO LeRa 
[eH BN = STN (OY + OK 


1954° [F OX < w THEN GN = PT + BN 
yaw OL = OF * BRel> * COS © ON) 

Lar’ Ge = OT + BRE Is = STN © ANS 

lesa OH = SHR 'O1 « 2 + Oe « 23 

18S Ol = SHR (OR «~ 2 + OF «~ 25 

time [IF OI > = TH THEN 1138 


l11iM O01 = OK:0e = OY 

ils BRET? = J + 1 

Lisa ‘YAC Ls = [NWT (Yee [lo + O24 
lisy secl2 = INT ¢8Be To * AL 
Lif NET 

L186 OR ETURN 
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Lira 
Live 
lirs 
i174 
Irs 
Lisw 
Lisw 
Ve ed 
Plas 
Lise 
11:35 
1136 
1197 
1138 
L2um 
i2ia 
1220 
12a 
L244 
{250 
eo t 
L2a2 
1 ew 
[254 
1255 
1255 
1257 
123s 
Se 
1238 
1235 
1234 
Se) 
13H 
i3lv 
et 1 
Pl. 
131-3 
i314 
$324 
oR. 
1344 
13548 
1 3au 
L378 


FEN 
REM 
BEM 
REM 
REM 
FAR I 
FOR J 
REM 
REN 
FEM 
FEM 
FEM 
FEM 
REN 
FOR I 
FOR J 


ATTRITION SUBROUTINE 


SERN GUT GTTRITTION 
WALWES FOR ALL LINTTS 


i Ta Merc l> = 4: 0RC I? 


= es NE 
i TO NeRFO J) = B30RC > 


4: NEXT 


ITHECK FRNGE BETWEEN 
IPPONENTS. IF RANISE 

<= 278 ATTRITION WILL 
NCCUR. REQUCE RATE OF 
RED MOVEMENT TO 2 44S 


WG) oy 


= ] 
= 1TO8N 


SBTC T,.J3 = B RTC I,[)>2 = 


PG = SOUR CC SBC TD — KREIS d «© Be PRE Tm PRE J od 
IF Ri5 % 2fbe THEN 1314 
PRC J) = 2 

REM 

REM INCREMENT THE NUMBER 

REM OF TSTS FIRED RT BY 

FEN EACH LINIT. COMPLITE 

REM RTTRITION CRUSED BY 

FEM FED ON BLUE ® WICE 

REM WERSRA 

FEM 
BRO [> = BFC IT? + 1 
RFej) = RFC J) + 1 

IF PROJ? = - & THEN BTV [,.7 = @: BOT 1295 
Bie l,J2> =~ Bw“ N © CL = (Po « SRC TI) © Re i) 
IF SReI> = 3 THEN PTC,[ 2 = e SATA 1318 
RT Je.DT> = BB“ M BCL = Cs «© SROI9S © SECTS 
MEST = NEXT 

REM 

FEN ‘SUM TOTAL STTRITION 

REM FOR ERICH LINIT 

FEM 

POR J = 1 TOM 

FOR J= 1 TO N 

IF RE“ > =H OR REC [> = & THEY ISSR 
MBCis = OCs + BE T,J7 « SRY 
DRC .J9 = OFC I> + RTC U,IT 3 © SPOT: 

NEST : MEST 
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bari oo REM 

Lara FEM COMPUTE NEM OUINTT 

1373 REM %‘STRENISTH. SUBTRACT 
1374 FEM UNIT ATTRITION FSM 
L375 REM CURRENT UNIT STRENGTH 
1376 FEM IF UNTT STRENSTH EOFS 
is’r 0 0OC REM) «6ONEDATIVE, MAKE [TFT 
13r°3. FEN EBUBL H. 

1373 REM 

1S3H FOR [ = 1 Te 
139M 56C1I> = SCT 
l¢mu 6o[F SCT) < 
i414 [F See[9 = = 
L38cn NEXT 

i430 FOR J=1 TO H 
[$a SRC I) = SRC.0 -— DR 1% 

1456 IF sRcJ> << & THEN SROs = 4 

346m [fF SRkeJo = = RM * Pit 1% THEM SPC 15 
147m NEXT 

1350 FETLIRN 

idsH FEM PRINTOUT -SLIMMOR'? 

1Smm IF PR = 1 THEN FPR# 1 

1914 PRINT Ti: POKE 36.7: PRINT Ree kk 13 
1515 POKE S6.ce:s PRINT RC KK 5 

1Scu PRR 

153@ RETURN 

1944 FEM %(‘SQUR FESULTS IN 8 FILE 

1554 PRINT 

1555 INPUT “On YOU WANT GO FILE OF THIS PINTS MOTH C'peN "SHS 
156m [FF AS = "N" THEN 1754 

157@ OS = CHARS (44 

i334 INPUT “NAME OF FILE FoR FORE | SHE) ORTOY “ars 

1834 PRINT OS; "OREM “sFS5". CORPS" 

16a «6PRINT OSs"WRITE “;FS5".C0EFS" 

15146 PRINT As PRINT B: PRINT Ps PRINT ff 

162 PRINT O3;s"CLOSE "GFS5". CORRS" 

16-3) SO PRINT OS;"0PEN “ZFS;". RESULTS” 

1634 PRINT O3;"WRITE “ZFS:". RESULTS" 

165% PRINT NK + 1: PRINT NK + 1 

lbeu FOR JT = 4 TO NK 

1b7H T =f G 

1686 PRINT Ts: PRINT Rie 7 

1BSH =~NERT [ 

lvoe FOR [ = 4 TH NK 

rig T=tT<3 
Lree PRINT Ts PRINT Rie [% 

Ps HET 7 

Lede PRINT OUSs°CLOSE “sPS5". FE SULTS" 
1P°S6 RETURY 


TA 
ie Hee [3 

HEN SR¢&TS = 

BM + Bile [> THEN BRC I 3 
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oe 
ree 
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FEM PRINTOUT-SUMMERY 
HOME =: WTRB S 

IF PR = 1 THEN PR# 1 
PRINT “TIME [35 “sT 
PRINT =: PRINT 

PRINT "“SLUE FORCES" 


PRINT “UNIT"s TABC 7 >; "¢CnORD" 


POR [ = | TOM 


PRINT " “sla TARC 75% 


NEMT 


PRINT : PRINT “RED FORCES" 
PRINT “UNTT"; TARE Fas"MCNORO": TOR? 


FOR [— = 1 TO WN 


PRINT " "sla TABC FISMRC Io; TABE 


NEXT 
PRINT + PRINT 
Pra 


INPUT “HIT RETURN TO CONTINUE" 


RETURN 


, 
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FORCE MAKER PROGRAMS 


i REM FORCE DORTA-FILE MHKEH 

24 0«OOIM BXC5,183,8'5,109>,R1° lo> 

34 6 6GOTO 176 

441 HIME 

Sa oINPUT “NeMeE OF FILE? “snes 

Ba sO =) OCHS OO 

f4y PRINT OS s"0PEN “SNSs"," SMS 

36 PRINT OS s"WRITE “GNSS. ard 

‘44 6 PRINT M 

144 FOR I = 1 TO M 

lift) PRINT Bicl? 

126 FOR J = 1 TO BiLT> 

13H PRINT 88CT,3: PRINT BY T.% 

134 NEXT J 

15H NEXT I 

ls4 PRINT OSs"CLUSE "“SNSs". "SME 

1rw HOME 

13a )06«6 TRB OS 

194 PRINT “MENUS: ' 

“04 PRINT "1 MAKE FILE FOR BLUE FORCES" 
214 PRINT "S MAKE FILE FOR RED FORCES" 
226 PRINT “4 WuIT" 

230 PRINT =: INPUT "WHICH? “32 

244 IF € = 4 TREN END 

254 IF © = 1 THEN M$ = "BLLIE" 

254 IF 2 = 2 THEN MS = "SEN" 

274 HONE 

234 PRINT “WHAT [5 THE FORE TSE OF “aise: TNPHIT "“ " GH 
23H FOR T = 1 70M 

set PRINT : PRINT “FOR UNTT “st 

319 PRINT “HOW MENY CONRRINATIMG POTMTS "5 
415 INPUT “S TNCLUOTNG STORT LOCATIONS oF "sR r4 


Reger = 1 Fis Bit 1. 


© a” 


436 FPeINT 

aah)06 PRINT “ENTER GRID COAPAS FOR CONPOINATING POINT = ° 5. 
224 INPUT “s-CnRO [5S “s6se Tot 

464 [NPUT “Yeninro [5 “sR Ti 

ay NERT J 

234 NEST I 
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LAAN os tn & : 
a te Oe, eS 
aSee Se Se Se oS oe Xe 


Ht 
«1S 


aX 


PRINT =: PRINT TREC t2o7"ee2 DATA CHECK see” 
PRIHT =: PRINT +: PRINT 

PRINT “UMIT PTs TREC Sas"8—nomPO"; TAR? IF o5°%SCnIIRo" 

FOR I = 1 TOW 

FOR J= 1 TO Beto 

PRINT TREC 23515 TARE Rigs TRAC 1G ssRRC TT, TRB Taoseee ld) 
NEST J 

MET | 

[MPUT “GN ISHONGES?T CeeN YG 3S 

IF 23 = “NM THEN oy 

PRINT 

INPUT “WHICH UNIT? ' st 

INPUT “WHICH POINT? “SJ 

INPUT “S=ZQ0RO IS “s88C 1.9 

INPUT “YeiDnnRO [5 “BY 1,139 

bHTOQ 3a" 
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C. RESULTS READER PROGRAM 


Los fri o--< 
1 a i. 


me nmy Tr ttt & & 
Se ase Se Coe Se 


oo * g-< | ened oye 
Cob Tale ee Es ee Ge 
ahr bye Xs 


&: 


Lis 
Lt 
43) 4) 


HME : THe 4 


OS = WCHRS (4% 


INPUT “HEME OIF FILE FOR FORCE |e fQTAy "ses 
PRINT OSs"GPEN “sFE&s" COEFS" 
PRINT OSs"FESO “sFss". Cners" 

INPUT Bs IMPUT 8: I[NFUT Ps: INPUT 8 
PRINT OSs "LUOSE “sFSs". COERS" 
PRINT OS;"0REN “ZF3Ss". RESULTS" 
PRINT OS;s"FREAO “sess". RESULTS" 
INFUT N: [MPUT M 
MK = N= 1 

DIM Bat NK D.RHC MK 

FOR I = 4 TO NK 

INPUT Ts IMFUT Bric To 

NEST I 

FOR T— = 4 Th NK 

INPUT Ts (HRT mee To 

NEST I 

PRINT O3s"CLOSE “sPs7". RESULTS" 
INPUT “PRINTER ONT “sHs 

MUME : UTTAR 3 

IF as = “Y" THEN FR 1 


PRINT “A = "cH 
PRINT “8 = "36 
PRINT "Po = “GP 
PRINT "Q = "GIR 


PRINT 
PRINT "TIME" cs: POKE GR.c: PRINT “SLUE": 
PUKE 2Rectts SPRINT "SED" 

FUR KK = @ THO NK 
T = SA = KE 

PMINT Ts: POKE =5,¢: PRIMT Bee KK 4; 

POKE 25,28: PRINT Reckk 2 

MexT KK 

PR i 

ENU 
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D. MAP MAKER PROGRAM 


REM HEP CATALCE MAKER.~FEADER 

MW DIM BLO ie, ly,le? 

3 HOME : UTAR 3 

44 PRINT “MENUS 

= PRINT “1 MAP MAKER" 

Ba OP RINT "2 HAP FERODER" 

‘QO PRINT “@& QUIT" 

34 PRINT : [HPUT “WHITH? “si 

34 TF # = @ THEN END 

Hy UN * BOTU L2b,.20H 

{1H FEM MAP CATALOG MRKER 

L246 HONE : WTRB 2 

tsa PRINT) TRBC 1133"*e* MOP MAKER +=" 

144 PRINT : FRIWT 

lS@ INPUT “NEW OR OLO MAP? CNY "Ses 

16h IF 235 = "N" THEN [I = trJJ = 1: GOTO Zen 

1°44 INPUT “NAME OF MAP? "sms 

Ls OS = CHRS 4% 

igH 6GOSLIB Ba 

“44 INPUT “ENTER STORTING SRID SPUGRE *Tl.ds “atta 
214 607T0 258 

229 INPUT “SIGE CF MAP ¢ leet HM OSDURRES OC "Glo 

234 INPUT “ENTER LOWER LEFT GRIOQ SGUARE CES. 56,359 “sleet 
236 [NePUT “BOSE ELEWATION OF MAP? " SHLO a. 

asa FOR ~ = [I TH L 

26H FOR J = JJ TO W 

S7MW PRINT : PRINT 

234 PRINT “FOR GRIO SQUBRE “sti"."5.J 

259 INPUT “NUMBER OF HILL MASSES? “SHLET.J.A KE = ALO TSS 
344 OTF KK 6= 4 THEN 2am 

218 PRINT “ENTER THE “GKEs" HILLS FOR SRID SAQUGRE “sla"." GJ 
se FOF & = 1 TO KE 

pew PRINT “RILE "sks" [3S “3: INPUT HLe 1.6 3 

336 HEAT = 

ao gj = f 


SoH PRINT : INPUT “SRE? 6 Yen os 
se TF ws = "Y" TREM shin 
23H HEST oJ 


43 
339 WERT £§ 
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Och & Gib 
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my Ted meg Wd oe 


f. US haw 
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ee 
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$6.8 


ROME 3: WTRB 4 

HPUT “NAME GF MAP FILET “ste 
OS = CHES ¢4> 

PRINT O3;"0PEM "sMS5". MAP" 
PRINT DSs"WRITE “sMSs". MRP” 
PRINT L: PRINT W 

PRIWT Liss PRINT L'? 

PRINT HLS 3.8.1? 

FOR [T= 1 TOL 

FOR J= 1 TO W 

PRINT HL [sJ.e? 

IF BLiI.wJ,9) = 8 THEN SSA 
POR K = 1 TO HLC I,J. 
PRINT HLilatsk 3 

MEST K 

NEST 

NEST I 

PRINT OSs "CLOSE “sM3i".MaP" 
PRINT : PRINT : PRINT : PRINT "... DONE": PSINT : PRINT 
INPUT “fi ‘YOU! WANT Th CONTINUE [NPUTTI MS KBTAY CHM’ “Ges 
TF F5 = "Y" THEN San 

IMTOO SH 
OS = THRS ¢4a5 

PRINT OSs"0PEN “shMSs". MGR" 

PRINT OS;"RERO “shMS;5".MHRP" 

[NPUT Ls INPUT W 

INPUT Liss INPUT LY 

INPUT HL a.,”> 

POR [ =1t TOL 

FOR 1 = 1 TOW 

INPUT HL T,J.8 2 

IF HULITsJ.4> = & THEN FSA 

POR K = 1 TO ALE I.Jok: 
INPUT HLS ToueK 


HET 
mest oJ 
NEXT I 
PRINT OSs"CLOSE "shS5". HAP" 
RETRY 


SEM MAP CATALOG FILE RESDER 

HIME : UTBE 2 

PEINT TREY tili1s"2e. HAP SEQNER «s«" 
PRINT : PRINT 
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ee va 
Ih On & Gi 
Se aye Se Se Se he 


ee Ce Wey We © ee 
fy: Ga = 


age Ub ig type ag Cpe og ve eye age 
i i> on & Gh es Le 
SADA DA DX AO DCU 


e- 

x 

Se 
aM] 


LBL 
Lian 
1.24 
Lad 
[Hse 
Laan 
Lar 
lyaw 
Lys 
Llu 
LLie 
Lice 


res 
cr, tye 
Se eS Se 


[ 
I 
Os 


PRINT OSs"0PEN "aiSa" 


p 
I 
I 
I 
es 


NPLIT “NeEME DF MAP FILET “sMs 
NPUT “PRINTER GN? CN “GAS 
HRS 645 


Map" 
RINT OSs"READ “shS5".NHP" 
NPUT Ls: INPUT WH 

MPUT Liss INPUT Li? 

NPUT HL 3.9.0» 

or I Lene 2 


FOR J= 1 TO W 


l 
: 


NPUT HL Ts ek? 
PF HLS Lod okt? 


H THEN : 


FOR K = 1 TH ALE Ts, 

INPUT Ale toJek 3 
NEXT K 
NEXT J 
HEXT I 

PRINT OSs"CLOSE "“shSs". MAP” 

HOME 

IF 4S = “Y" THEN PR 1 

FRINT “NAME OF MAP: “3MS 

PRINT "“LUWER LEPT ISRID SHUBRE: “sles” 
PRINT "SIZE OF HRP: "sls" BY “GW 
PRINT “SASE ELEVATION: "SHLO AAAS 
PRINT 

FOR IT =t TOL 

FOR J= 1 TO W 

PRINT [.JeseUt I,J, 


FOR K = t TO BLOT... 
PRINT HLilsJow oa" 8G 
NEAT K 
PRINT 
NEXT J 
RT 
PR 1 
Ren $s ee iNT 

(NPUT “HIT RETURM TO CONTINUE” 25 


BOT 3h 


HEXT [ 
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HILL MAKER PROGRAM 


ee tt MILL OHTS FILE MAKER-RESMER 
24 GIh Ae = ee Jee Cis 3 


a6 Ki 12 = SskCS5 = SekK£ 33 = toekC do = Te 
$9 KOS) = S2rKCR) = 24eki ro = 29s Ki) = 32 
34 REM HILL GHTR FILE N@KER 


bd HOME : UTAB 2 

ra PRINT “MENLs 

34 4 PRIHT "1 HILL MAKER” 

34 PRINT “2 HILL RESOER" 

14H PRINT “8 BUIT" 

{la PRINT =: INPUT “WHICH? aH 

126 IF * = 8 THEN ENDO 

13H ON & BUTO 1358,3826 

144) FEM HILL DATR FILE MAKER 

1S@ HOME : UTRB 2 

16H PRINT TREC 1135" See HILL MAKER see" 
tra PRINT : FRIWT 

Ga OT MPUT “NEW OR COLO FILE? OMe “pee 
136 IF ad = "N” THEN 234 

264 GHOUSUB SAB 

S18 INPUT “WHICH HILL Ot YO) WeaNT TO STRRT WITH “aT 
224 I3HiTO 2650 

236 HOME =: UTHR 4 

2446 i] = i 

2S (UMPUT “NUMBER CIF HILLS?’ “aN 

cBY FUR T = II TQ H 

27k PRINT "ENTER THE OATS FOR HILL “al 
2s INPUT "SEOURD [5 “sHo l,l: 

“30 IMPHOT “PiOcRO IS “sHOT.2 3 

3446 INPUT “PEGK ELEWATION [SS “see [.3% 
sig INPUT "GRIENTATION ANSLE IS “sHET.4% 
328 INPUT “ECCENTRICITY Is “seHOILS) 
.c40 6 INPUT 6“SPRERD [5 “SHO 1.52 

344 LNPUT “Hex HEIGHT [5 “see I. 

204 INPUT “CUT HEIGHT IS “sHe [.32 

aba PRIHT =: PRINT : INPUT “SRUE? CYeNo "ses 
avA LF 23 = “Y" THEN =e 

254 NET I 

230 HOME = UITAR 4 

44h OS = (CHR 64% 

410 INPUT “NBME CIF FILE? "ans 

426 PRINT O3;"0FPEM “sNS3" HILLS, Lid" 
$30 PRINT O33" WRITE “sNSi". HILLS, Fe" 
46 PRINT 4 

$54 FOR T = 1 TON 

46H PRINT OSs"NRITE “ENnSs" HILLS, ost 


V2 





mine in & + & 
WR om) ich Ri Go On 4 
ra Bagg 


at 
Sy Ge Se Ss Ol 


ao 


rx 


os a] 


Nn 
‘Ts 


o™] 


on OS 


rin & 


oye af) se Seen ; 
Sets 
Sai ie 


2B 
Ge Sp a 


FOR J= 170 & 

PRINT HET.) 

HET 

NEST I 

PRINT OSs"CLOSE “sNSs". HILLS" 
PRINT : PRINT : PRINT "... QGUNE" 
PRINT : PRINT 

INPUT "Co ‘PON WENT TO CONTINUE TNPLITTING DOTA? Cle No SAS 
IF eS = “Y" THEN 214 

ITO SY 

HIME 
OS = CHAS ¢ 4) 

INPUT “NAME OF FILE? “sNs 

PRINT OSs"CIPEN “sSNSs" HILLS. Lae" 
PRINT O$;"REQD "SNSs". HILLS. Fa” 
INPUT N 

FOR I = 1 TQ N 

PRINT OSs"FESBO “sNSs". HILLS. Rat 
PUR J =1 TH 8 

INPUT HE TU? 

NEXT J 

MEST I 

PRINT OS3s"CLOSE “sNSs". HILLS" 
KUNE 

PRINT “HUMBER OF HILLiSs “sNe PRINT 
Ma = | 

FOR I = 1 TO N 

IF M = 11 THEN PRINT :M = } 
PRINT I3 

FOR J=1 TQ & 

PRINT The KC Josue TS t3 

NEXT .J 

KRINT 
MoM + I 

NEXT I 

RETLURMM 
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32 
sh 
a5 
36a 


ara OS = 


334 


yay 
LylW 
Wee 
1W38 
a) 
1yWS8 
1 Yee 
Lar 

Lys 
1ygH 
Liga 
Lite 
Liew 


Se eye oy 
Me iS iS 


REM HILL ORTA FILE FESTER 
33H rR = 4 
MEME : UTAB = 


PRINT TABS 


PRINT : 


PRINT 


THES 6 $5 


{A234 HILL READER see" 


INPUT "NAME OF MAP? “GN 
INPLIT “PRINTER ON? C/N) “52s 


IF £3 = 


"y" THEN FR = 1 


PRINT OSs"0PEN "GNSs".HILLS,» 134" 
PRINT OS;"RERO “GNSs".HILLS. RO" 


INPUT N 
FOR I = 


1 Tf) N 


PRINT OSs"READ "SNSs". HILLS. RUGT 


FOR J=i TH & 


INPUT HOT oJ) 


NEST J 

NEST I 
PRINT 
HUME 
lear 
PRINT 
PRINT 
PRINT 


KK = kta + 
IF PR = } 
PRIHT TRBe k 


NEST 
PRINT 


UerelLUSse "sees. HELLS” 


e 
¢ 


1 THEN PR# 1 
PRINT NSa" HILL DATA": PRINT 
"MUMBER QF HILLS: “GN 


2 


cy) aa le a 


NEXT I 
PRe A 
INPUT 
SOTO 5 


« 


EM PRINT <M = 1 


ad 
THEN POKE GA.KK:s PRINT HOT. cs2 BOTA Lise 


Soule Maret pel oles 


"HIT RETURN TO CNTTNUE” i5 


val 
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F. 


WOODS MAKER PROGRAM 


REM 


Gi FLY lie LY eS | 


HOME 
PRINT 


: WTAB S 
“MENU: " 


PRINT "1 HUODS MAKER" 


PRINT 
PRINT 
PRINT ; 


[NPT 


“WHICH? 


[F * = 8 THEN END 
BN « SOTO Len. rr 
REM WOODS CATALCIS MAKER 


HOME : WTHB 2 

PRINT TAC t13s"se2 WOODS MAKER ++" 
PRINT : PRINT 

INPUT “NEW OR OLO MOP? ENED “sss 

IF 2s = Mw hytt THEN ag = i ane = ts STE 224 


"2 WOWDS RERCER" 
“9 MUTT" 


is * t,t 
aa 


TNPUT “NAME OF MAP? "GMS 


O63 = 
BLISUE Bett 


HRS 643 


INPUT “ENTER STARTING SRID 


SOT 24a 


HQ003 CATALOG NAKER-RERDER 


SMQUIGRE ¢T -.l3 


[HPUT 
INPUT 
FOR I 
FoR oJ 


PRINT ; 


PRINT 


"SIZE OF MAP « tame 
"ENTER LOWER LEFT 
Ll ee 

JJ TO W 

PRINT 

"COR GRIO ScUARE ": 


M 
3k 


SMUBRES °? "aL 


ID SQUARE “EG. SH.53: 


o 3 if . | 
4 # Pm 


IHPUT “NUMBER GF FORESTS? "GFLY[.J.4%2kK = 

IF KK = 4 THEN [YA 

PRINT “ENTER THE “3kKK5" FORESTS SNR GRIT SaLlor 
FOR Kk = 1 To KE 

PRINT “POREST "sks" [3S “as: INPUT FPL Tsu.k % 
HEMT K 
vo | 


PRINT : 
IF G3 “ bey te 
NEXT J 
Meal [ 


THEN S48 


BRS) 


INFUT (SSGQUIE? 6 pup oo" sos 


: 3 if I ao 


m 


a oe 
er - 8 


ou 





“OA oH ON 
fra oe Re 
Ss Ss 


my A i 
Gh: 


Io ON 


oe oS CS 


a cy Te a Ts Ts 8 Ss Te iy td CH 
a 25 0 eit & Life = © co Os ~. 1s 
0 oS Se ce cy Se Ee Ge Si Se Ee CS Se eS 


2 teh Fed es 
61% c& 


iN f tel 
Se Si Se 


=) =) ~ J am | =a ~~) oF 


HOME : VTHR 4 
INPUT “NAME DF MAP FILE? "GMs 


aus = CHRS (43 


PRINT O33"0PEN “sMSs". WO00S" 
PRINT OSs"°WRITE “shSs".WO005" 
FRINT Ls PRINT W 

PRINT Lis PRINT L'? 

FOR T= 1 TOL 

FOR J= 1 TOW 

PRINT FLY I,J. 2 

IF PLil.J,.4> = @ THEN 358 
FOR K = 1 TH PLO TsJet2 

PRINT FLO T.J.Kk 3 

NEST K 

NEST OJ 

NEXT I 

PRINT O3s"CLOSE “aMss". WOODS” 


PRINT =: PRINT : PRINT : PRINT "... 
[INPUT “OO ‘OL! KENT TD CONTINUE THRPHUTTIHG GRETA? Cleon od MGS 


IF 23 = "9" THEN San 

ut 3a 
OS = CHRS “42 

PRINT OSsS"0PEN “GMS; ". WOnDS" 
PRINT O$;"RERD "GMS". NO00S" 
INPUT Ls INPUT W 

INPUT Liss INPUT L'? 

FOR [T = 1 TOL 

FOR t= 1 709 W 

INPUT FLO T.J.2 

IF FLil.J,9°> = & THEN Fey 
FOR K = 1 TD PL T pd stt 

IMFWT Fluo ledek 3 

NEST K 

NEMT J 

MET 

PRINT OSs"CLUSE YM sMsa". HOODS" 
RETLIRM 
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PRINT : 


PRINT 
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REM HUUDS CATAL FILE FRESDER 


HOME : WTRE 2 
PRINT TRBC Lijs*ee*e WOODS REGQER +.+" 
PRINT : PRIWT 
[INPUT “NAME CF MBP FILE? “sms 
IMPUT “PRINTER ONT 6 YON “GAs 
CHRS (32 
N33 "CiPEN 
0$5,"RERQ 


PRINT 

PRINT 
INPUT Ls: 
INPUT Lis 

FOR IT = 1 

FOR J= 1 TOW 

INPUT Fle Il oJ se % 

IF FLO I.J.4)9 = 4 THEN 3an 

FOR K = 1 TIX PLOT. 

INPUT PLel,J.K3 

NEWT kK 

NEAT J 

NEwT f 

PRINT OSs"CLOSE “iamMss".W000S" 

HIME 

IF 8S = “Y" THEN PRE 1 

PRINT “NEME WF MAP: “5M 

PRINT “LOWER LEFT GRID SuRRE: 
FRINT “SIGE WF HAPs: 
PRINT 

POR [T= 1 TOL 

FOR J=1 TOW 
PRINT [,J.Fit [.J.4? 
IF FPLC T.J,e> = 4 THEN 111h 
FOR K = 1 TY FLO T.J.4: 
PRINT PLC IT sdeK os" “a 

MET K 
PRINT 
NEST J 
PRINT : 
PR# 
PRINT : PRINT 
INPUT "HIT FETURN TO CONTINUE” iY 
Bit 2H 


“SMss" WROS" 
“sMs;" WOODS" 
INPIT +4 

INPUT Li? 

wee 


“sls” BY “GH 


NEXT f 


re) 2) 


Bates ly 





G. FOREST MAKER PROGRAM 


ltt FEN FOREST MAKER-PRERTER 
24 OM Fe leb,BI,KiR 2? 
49 KCl) = SskC29 = leskt > = Deki ad = 24 

$4 K05> = TO:k¢ Bho = GB 

3 ANME : WTAR 2 

mY =60PRIWT “MENU: 

‘4 PRINT “1 FOREST MAKER 

34) PRINT “2 FOREST RERDER" 

38 PRINT "@ QUIT" 

144 PRINT : INPUT "WHICH? “Ge 

{ig IF * = & THEN END 

126 ON & GOTO 144,20 

13H «OREM FOREST OQATA FILE MAKER 

is HOME : LWTGB = 

[Sa PRINT) 60TARE Litas" see ENREST MORKER +e<" 
164 PRIHT : PRINT 

Lr INPUT “NEW OR OLO FILE? “M-N> "GAs 

13H #xIF 83 = "N" THEN 2246 

igy BOSUB Saa 

266 INPUT “WHICH FOREST Of ‘POL! WANT TT START WITH? "ETT 
214 mHogTd 259 

ecu HOME : VTAR 4 

a3 [I = 1 

234 INPUT “NUMBER OF FORESTS? "GN 

254 FOR I = II TON 

254 PRINT “ENTER THE ORTR FOR FOREST "37 

sv [MPUT “SCCIRD [5 “see l.i3 

234 INPUT “YCrmRD IS “Zee l.2% 

234 INPUT "HEIGHT OF TREES [5 “sPel,3% 

240 [NPUT "“SRIENTATION SNGLE IS “cFé..4% 

Sig INPUT “SEMI-MQJ0R BKI'S LEMSTH Ts “Gee l.a% 
se INPUT “SEMI-MINCR AXIS LENMSTH [5 “sFe Ta: 
339 PRINT : PRINT : INPUT “SRUE?T CYens "Zest 
434 IF 235 = "" THEN An 

356 NEXT I 

4268 LIME ; UTR 4 

374 OS = CHRS (4% 
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[HPUT “MAME OP FILE? “ans 
PRINT OS;3°0PEN “sngs" FORESTS. Las 
PRINT OS3;" WRITE “sNS35",FORESTS,. BH” 

PRINT WN 

FUR IT = 1 TOW 

PRINT O3s"°HRITE "“sNss" FORESTS. R'iT 

POR | = 1 TO 

PRINT Fe Tso 

HEAT 

Mest) | 

PRINT OSs"CLOSE “iands". FORESTS" 

PRINT : PRINT : PRINT “... GONE" 

PRINT : PRINT 

INPUT "00 YOU WANT TO CONTINUE TNPUTTING MATA? C'YeN Go “GAS 
IF 43 = "YY" THEM San 

iKsiTO Se 

ACME 
OS = CHES (4; 

INPUT “NAME OF FILE?’ “SNS 

PRINT O$2"QPEN “GNSs" FORESTS, 1.3” 
PRINT GS;s"REHD “sNSs". FORESTS, Fn" 
TNPUT N 

FOR I = 1 TON 

PRINT OSs"RERD "ZNSs". FORESTS, FR" :I 
POR $= 1 7TO85 

INPUT Fed. 


NEXT J 

NEXT [ 

PRINT O8s"CLOSE "“iNSs". FORESTS" 
HOME 


PRINT NBs" FORESTS LIST": PRINT 
PRINT “NUMBER DF FORESTS: “an: PRINT 
oO pee 

FOR I =i TON 

IF M = 11 THEN FRINT = = } 

PRINT [3 

FOR § = 1 TOS 
PRINT TARE Ke Jo sFe Tn 53 
NEST 
PRINT 
Mat + ft 
NET I 
PETURM 
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REM FOREST ORTH FILE FPEADER 
HOME : UTHB = 


320 PR =o 


PRINT : PRINT 


35a 05 = CHRS (42 


"; og 
& 


6 CEs oe Ct ee 
oS oa Se SH 


ro os St 


INPUT “MRBME OF MAP? "SNE 

INPUT “PRINTER ON? (AND "ses 

IF 23 = “Y" THEN PR = 1 

PRINT O3s"0PEN "sNSs". FORESTS. Lag" 
PRINT OSs"REAO “sNSs". FORESTS. RQ" 
INPUT N 

FOR [ = 1 Ti N 

PRINT O3s"READ “SNSs". FORESTS. RP’ Gt 
POR J = 1 TO 5 

Rae Fs ls 2 


MEXT J 

NEST I 

PRINT OSs"CLOSE “sNS3". FORESTS" 
HIME 


[IF PR = { THEN PR# t 


PRINT : PRINT Ns" FOREST ODRTA": PRINT 
PRINT “NUMBER UF FORESTS "3N 

FRINT 

M=1 


FOR IT = 1 TO MN 

IF 4 = 11 THEN PRINT :M = 1 

PRINT I; 

FOR J=1 TO 8 
KK = kt.J>o + 3 # J 

IF PR = 1 THEN POKE 36.KK: PRINT Fetus: GATTO t1l2e 
PRINT TREC Ke Jsa5FOT ,33 

NEST J 

PRINT 


She > 1 


NENT I 

Pes tt 

TNPLUT “HIT RETURN TO CONTINUE” is 
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1 FEM UPOATED 2 FEB Se 

“4 FEM OGAMSTAC TERRAIN MHOCEL 
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OIM 2BeS:,2R¢5°,HL 14,106,153 | 

$4 OIM BPCS>,RPCS 1.5865 2,.5R65 0,29 5,14% 

249 OIM BFOS7,RFC59,67T(5,52,8T65,5 29,0863 >.,9R¢ 5 3 

BY GIM BRC 5S),RRCS2 

‘4 OIM B1¢52,R1¢5°% 

36 T = 4:0T = SaenT = 150 

= LHRS (4% 

P= W390 = 4.o3 

114 PI = 3.1415926554:PR = 

Hz 2.3:38 = AWS 

W.S32:PE = 4.3 

~O2RM =, 
mi) 
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165 22 = 1 
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el 
ci 
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_A Qa 

Sf tr’ a 
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tf UIM Batt NK .RACNK » 

171 HOME 

Lr2 INPUT “NeEME TF MOP FILE? “ss 
ifs PRINT O3;,"0PEN “;MS;5".MaAP" 
L°+ PRINT O33"SEAD “shs5".MAaP* 
t1?S INPUT I: INPUT wW 

ivS [NPUT Lit: INPUT 1 

ter INPUT HELD 4,4,0% 

leis FOR IT = 1 TOL 

1°32 FOR J =1 TD W 

lsH INPUT HLO 1,0, 


L3L IF RLOL,J,4°5 = & THEN Lied 

82 FOR K = 1 Th HL TJ. 

LSS INPUT HELE I,J.k 3 

lo34 NEXT : NEIT = NEXT 

135 PRINT OSs"CLUSE “sMS5“.MHP" 

136 LS = LS - ltl? = Lit = 1 

Lisf PRINT O35"0PEN “SMSs".HILLS. Iso" 

138 PRINT OSs"RERO “snss" HILLS. Ra" 

133) TNPUT Nee OIM He NH, 5 

Su FOR [ = 1 Tl NH 

i131 06 PRINT OF;s"REQO “SMSs" HILLS. YG 
ig2 POR J = 1 Th & 
145 [HPUT He T..J 3 

(194 NEXT : NEXT 
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PRINT O3;"0FEN “GNSs". BLUE" 

PRINT O3;"RE80 "“sN3;", BLUE" 
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FOR I = 1 TOM 

INPUT Bl 

FOR J= 1 70 Bl 

INPUT BKC I,J3: INPUT BY I,J2 

MEST 3s NEXT 

PRINT OSs"CLUSE “sNSs". BLUE 

FOR I = 1 TO MeX8BCT® = BkKCLosL eB ls = Bre T,tasBee ta = 1 
PRINT 81: FOR JJ = 1 TO Bt: PRINT BKC Toto Toe NEXT 
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INI5UB Ls6B 
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PRINT "ELEVATION OF BLUE "sls" [5S "32 
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NEXT 
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PRINT OS3;°0PEN “GNSs"*. RED" 
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FOR I. = 1 TOWN 

INPUT Fi 

FOR J = 1 Th Ri 

INFUIT RECT. c: [INPUT Rod T,.52 

MEAT =: NEXT 

PRINT US;s"CLOSE “GHS;". RED" 

POR [T = 1 TO Ne#8ROTo = SXC TT.) o2'PRO To = BY l,i she To = } 
PRINT “SIZE OF PED UNIT “sla: TMRLIT SRe rs 


Picts = SRCT > 
RRIC HM = RACAD + SR I % 
ARE I > = 5 


SRC Las = ‘vee To 


ZR a = TNT «eRe lo - [eta = Liked = TNT Cpe To Lees = LY 


SUISUEB [Sey 
ee( la se 

PRINT "ELEVATION OF FEO “ZIs" [5 “ae 

POR [1 = 1 TQ HLE Ma, Wa, 92H? = HL kd od, Tt bene He. = ts NENT 
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238 PRINT = INPUT “fi ‘ool WENT THe SRTNTEO Cie een GS 
535 IF aS = "''" THEM FR = 1 

344 HUME : LTHe 3 

354 IF PR = 1 THEN PR# t 


364 PRINT 

SV) PRINT “TINE“s: PORE SA.72 PRINT "SLE" ss OnkE SR,2m%° PRINT “RED” 
Sau PR# 

739 IF 22 = 1 THEN %j(.BOSUB 17°6@: SOTO 81m 

RM 650SUB 1430 

Biy FEM MAIN PROGRERM 

bew T = T + OT 

o34 [F T < = MT THEN 54h 

335 PRINT “ENO GF BATTLE OWE TH TIME": 0SHIR 154@- ENT 
Bay =650SU6 734 

how G6USUB 1178 

BBY KK = KK + |. 

Bra BYCKKD = AsReC KK = 4 

boy FUR IT = | ‘TO MM 

“Ja BHC KK? = BUC KK + SBE I 3 


raw NEXT I 
‘lw FOR J = 1 TON 
fey RACKK > = RACKKD + SRO 


foM NEXT J 

oy IF 22 = 1 THEN GOISUB 17Sae GTN TRA 
Yo” iEnOSuB 1435H 
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Bo PRINT : PRINT “ENO OF BATTLE GUE TO STTRITION OF SLUE" 
hs %'sUiSHWB 134: END 

re IF RMACKK: > RE + Rada THEN Fea 

reo PRINT =: PRINT “ENG CF RATTLE MIE TO ATTRITION OF RED" 
res GUSUB 1544: END 

ray HOT Rew 

P34 FEM MNVEMENT SUBFROLUTINEG 

syy FOR I = 1 TOWN 


s1H J = BR To 

s20 OY = RY. + 12 = PRE T 

338 08% = REC L,I - La —- KRY Ts 

Ou LF Live = M4) MNO Cty * Mt THEN Gh) eT a “Je ea ty 


sa IF 0% = 8 SND OY < «@ THEN ON a ¢*¢ FT - 2: SOTO Sey 
6h HN = net ‘iy “ Fis 


sya [FF OX < W THEN ON = PIT + GM 
sa4 01 = aT * WRC I + mrs Yop 

334 U2 = OT + RRCIT? © SIN © BN% 

WW4 0H = SOR (O01 « 2 + fe « 2% 

Sig Ol = SOR (OK «~« 2 + OY « 3% 

teu IF OF > = GH THEN 358 

330 G1 = 08:02 = OY 

34M FPCT Oo = J + } 


163 





On i) 
fi 


e 8 
. 
SS Sa 


eg. coe Coe ce Ue 
‘ Fi gy a. ) 
Mon gf GIihie 


ge tg ogo 


6 ei a ie 
“et OR me Ss ThA 


m o- ifs 


eo Br [ 
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NEST 

FOR I = 1 TQ 


[ 
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mH ri 


' Sa se ‘a a 


Pl “ 2: BNTO 198a 
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D1 = OX:02 = OY 

vee 
Foi (Selo + oS) 
INT ¢xXB¢7T) + O19 
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NEST 
RETURN 
FEM 
FOR | 
FOR .) 
FHF J 
FOR J! 


a a) 


TTRITION SUBROUTINE 
1 TO MsBFO Is = a:0BCI2 = &: NEXT 
1 TO Ne REY J) = Ae ORe > = He NEXT 
1704 
Peeters 


24 BTi l,j) = @:RTCJ,0>9 = 4 


AAG = SHR COMBE TY — ERO Idd 
IF Ris > 27a THEN 1214 
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A RRC TO = 2 
MEM (CHECK LINE CF SIGHT 
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me = INT C#RO Jo © Laem = Lite eo = TMT oR da & Lee =K LY 
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L2he GBOSUB 1364 

1263 NEST = NEXT s 

1268S FOR T1 = 1 TH NHsHiTi.a@> = As NEXT - 
1265 FEM Lis [5 Ok 

l2r4 BROT) = BFC TD + 1 

L234 RFC J) = RPO do + Lt 

2S BTC1.7 ts W #¥ C1 ~ CRO A SEC TOD e SRE) 
1369 FTC JsT3 = Be fl = fo « SRO) * SRC TO 
1318 NEST =: NEXT 

is26 FOR I = 1 TOM 

igsH FOR J = 1 TON 
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1374 NEST : NET 
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195% IF SFeJ: < 8 THEN SRC Jo = 4 

Idem IF SFkeIo < = RM * Sie 1) THEN BRET) = = & 


1974 NENT 

14s KETURN 

13994 FEM FPRINTOUT~SUMMARY 

1Sy4 IF PR = 1 THEN PR# 1 

15Siy PRINT Ts: POKE 36.7: PRINT S&C KK 35: POKE SR.2&: PRINT BMC kK 
1528 PR 4 

1530 RETURN 

153s REM SAVE FESHLTS IN A FILE 

1354 FRINT =: [MPUT “OO YOU WANT RA FILE OF THIS RUN‘’S DATRY & eon 
1564 [F 435 = “H" THEN 17°54 

lof OS = CHRS «4% 

15s INPUT "NAME TF FILE FOR FORCE LEVEL DATA? “srs 
S36 «PRINT O3;5"0RPEN "“SFS;5" NERS" 

lenty PRINT OS;s"WRITE “;F35".C0ERPS" 

161i PRINT As PRINT 8: PRINT FP: PRINT 

ibe PRINT USs"CLUSE “sFS5".COEFS" 

16-26 PRINT OSs"0REN “ZFSs",PESULTS" 

le34 PRINT OSs"WRITE “GF35". RESULTS" 

i6bSH FRIWNT NK + i: PRINT HK + 1 

lobby FUR [ = 4 TO NK 

IBY Tse 3 

Lesh PRINT Ts PRINT Sec [3 
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FOR I = 4 Tl NK 

sel = 4 

PRINT Ts PRINT Reet 

MENT i 

PRINT O3s"“CLOSE “;F$3".RESULTS" 

RETLRN 

REM PRINTOUT <“SUMMARY 

HOME =: UTRB 3 

IF PR = 1 THEN PR ! 

PRINT “TIME [5 "“s3T 

PRINT : PRINT 

PRINT "BLUE FORRES" 

PRINT "“UNIT’ 3 TABC Pos"¥CONRD"; TARC 159."Y"YCRORO"; TAB & 
FOR I = 1 TQ 

PRINT “ “sis TREC P3OsKBC IDs TABS 15 3s¥8C 055 TARC 24 935R¢ I. 
NENT 

PRINT : PRINT “RED FORCES” 

PRINT “WUNIT"; TARC 72s"°¥COORO0";: TARC 150; PCONRE" > 

FOR [— = 1 TO HW 

PRINT “ “Sis THC PosXROT 3s TOBE 15 05'°RO To; TORS 24 5955R¢ 19 
NEST 

PRINT =: PRINT 

PR 

INPUT “MIT RETURN TO CONTINUE" :2s 

RETURN 

REM ELENATION COMPLITAT ICN 
CR = = 99999999 

= ALO YA 


FOR Ht = 1 TO HLS Xd.¥4,03 
He = BLO X4,'2?4,H1 2 
IF HCH2,.42 = 1 THEN 2156 


mS = % ~ HOCH2,1) © Lies = 4 — HOH? ,.2 * LM 
HL = Lom CHCHS,79 + eee? ~ SH )3 
Rl = Al # ¢HOH2,599 « 2 
Mil = HoH? 40 # PT ~ Lan 
a = HOH2,59 «~ 2 
arose (el =! COS (Cla a 2 + Bl eo SIN CCL) a BY Ce 
ae eee i Ce? FRI Se OA Tl A 2 le 
MS =co © COS (Cid * SIN (Ot? * (RL = 819% - Te 
[F HH2.°3 > HOH2.390 THEN CR = LOG CC HOHE. PD - AO He 
Wl = Pl *# #35 « 2 + PO & PS A EB + PR Se YS & YS 
He Me phe = i 
IF WI < ie THEN 215% 
PL = HHS, 3" + HOHS,F > 4 0 EXP COL a = 13 
IF Fl > 2 THEN = = FT 
MEET 
RETURN 
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